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34 PUOFESSOtt W. M. HICKS ON VORTEX MOTION. 

The chief part of the following investigation (Sects, i. and iii.) was undertaken with 
the view of discovering whether it was possible to imagine a kind of vortex motion 
which would impress a gyrostatic quality which the forms of vortex aggregates 
hitherto known do not possess. The other part (Sect, ii.) deals with the non~ 
gyrostatic vortex aggregates, the discovery of which we owe to Hill,"'" and investi- 
gates the conditions under which two or more aggregates may be combined into one. 
It is shown that it is allowable to suppose one or more concentric shells of vortex 
aggregates to be applied over a central spherical nucleus, subject to one relation 
between the radii and the vorticities. In all cases the vorticities must be in opposite 
directions in alternate shells. The special case when the aggregates are built up 
of the same vortical matter is considered, and the magnitudes of the radii and 
the positions of the equatorial axes determined. The cases of motion in a rigid 
spheroidal shell and of dyad spheroidal aggregates are also considered. 

The chief part of the paper refers to gyrostatic aggregates. The investigation 
has brought to light an entirely new system of spiral vortices. The general con- 
ditions for the existence of such systems, wdien the motion is symmetrical about 
an axis, are determined in Sect, i., and are worked out in more detail for a particular 
case of spherical aggregate in Sect. iii. It is found that the motion in meridian 
planes is determined from a certain function \jj in the usual manner. The velocity 
along a parallel of latitude is given by v ~f($)/p where p is the distance of the 
point from the axis. The function t//, however, does not depend on the differential 
equation of the ordinary non-spiral type, but is a solution of the equation 



d-ylr , 1 cPylr cot 6 d4r .,_ r df 

^ r 1 d0* r l dd H J df 



- •> 



where F and f are both functions of i//. The case F and fdfjdf both uniform is 
briefly treated. It refers to a spiral aggregate with a central solid nucleus, and 
is not of great interest. The case F uniform and fee xjj is treated more fully. If 
' f = hjf/a where a is the radius of the aggregate 

The most striking and remarkable fact brought out is that with increasing para- 
meter X. we get a periodic system of families of aggregates. The members of each 
family differ from one another in the number of layers and equatorial axes they 
possess. I have ventured to call them singlets, doublets, triplets, &c, in contra- 
distinction to the more or less fortuitous and arbitrary compounds dealt with later, 
and which I have named monads, dyads, triads, &c. Of these families two are 
investigated more in detail than the others. In one family (the X 2 family) all the 
members remain at rest in the surrounding fluid, In the other (the X 2 family) the 

* u On a Spherical Vortex," s Pkil. Trail a./ A, vol, 185, 1894.; 
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distinguishing feature common to all the members is that the stream lines and 
the vortex lines are coincident. 

The parameter X defines the total angular pitch of the stream lines, on the outer 
current-sheet, viz., up the polar axis and down the outside ; although in the aggregates 
with more than one axis these lines are not one continuous stream line. The first 
aggregates — with \< 5*7037 (the first X 2 parameter) — behave abnormally. Beyond 
these we get successive series, in one set of which the velocity of translation is in the 
same direction as the polar motion of the central nucleus, in the alternate set the 
velocity is opposite, and the aggregate regredes in the fluid as compared with its 
central aggregate (see fig. 3, Plate 1). The physical analogue of these aggregates 
is obvious. It is specially enlarged upon in the abstract.* 

Suppose we set ourselves the problem of making a set of aggregates with greater 
and greater angular pitch. As we do so w T e shall find that as the pitch increases 
the equatorial axis contracts, and the surface velocity diminishes. On the outer 
layers (ring shaped) the spiral is chiefly produced on the inner side facing the polar 
axis, until on the boundary itself the stream lines flow in meridians, and the 
twist is altogether on the polar axis. The pitch can be increased up to a certain 
degree. As this is done, the stream lines and vortex lines fold up towards one 
another, coincide at a certain pitch, and exchange sides. When an external angular 
pitch of about 330° is attained it is impossible to go further if a simple aggregate is 
desired. If a higher pitch is desired it is attained by taking it in two parts. First, 
a central spherical nucleus of the same nature as the former, in which a portion of 
the twist is produced, and outside this a spherical shell, in which the spirals have 
the same direction of twist, and complete the pitch to the desired amount but in 
which the spirals are traversed in the opposite direction. With increasing pitch 
this layer becomes thicker, and its equatorial axis contracts relatively to the mid- 
point of the shell until another limit is reached ; the stream, and vortex lines again 
fold together, cross, and expand as this second limit is reached. If a larger pitch still 
is desired there must be a third layer, and so on. The first coincidence of vortex 
and stream lines takes place for an aggregate whose pitch is 257°'27 / . Whenever a 
maximum pitch is attained the aggregate is at rest in the fluid. This is first 
attained for an external pitch of 330°'14'. Beyond this there are two equatorial 
axes. For an external pitch of 442 *37' the stream and vortex lines again coincide, 
the internal nucleus gives 257°* 27' of the pitch and the outer shell the remainder, 
and so on. 

At the end a theory of compound aggregates is developed similar to that in 
Sect, ii for non-gyrostatic vortices. It is not worked out in detail in the present 
communication, but the conditions are determined for dyad compounds, whilst a 
similar theory holds for triad and higher ones. Each element of a poly-ad may consist 
of singlets, doublets, &c. The equations of condition leave three quantities arbitrary— 

* ' Roy. Soc. Proc.,' vol. 62, p. 332. 
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as, for instance, ratio of volumes, ratio of primary cyclic constants, ratio of secondary 
cyclic constants. The full development of this theory is. however, left for a future 
communication. It is clear that spiral or gyros tatic vortex aggregates are not confined 
to forms symmetrical about an axis. Their theory is however much more complicated. 
If we take any particular spherical aggregate with given X and primary cyclic 
constant (/x) s the energy is determinate. We may, however, alter the energy. If it 
be increased, the spherical form begins to open out into a ring form, whose shape and 
properties have not yet been investigated. If the energy be increased sufficiently 
the aperture becomes large compared with the thickness of the rotational core, and 
approximate calculation can be applied. The differential equation for xfj is given in 
Sect, i., but its development is left for a future occasion. After that I hope to deal 
with the question of stability, and then more fully with that of the conditions of 
combination. The new field opens up so many questions of interest that other 
workers in it are welcomed. 

Section i. — General Theorems, 

1, To give an idea of the nature of the motions considered in the present investi- 
gation, consider the case of motion of an infinitely long cylindrical vortex of sectional 
radius a. The velocity perpendicular to the axis inside the vortex will be of the 
form v =f(r) 9 where f(0) = 0. Outside it will be given by v = Va/r, where Y=f(a). 

We may, however, have a motion in which the fluid moves parallel to the axis 
inside the cylinder with rest outside. The velocity will be of the form u = F (r) 
inside, where F(a) = 0, and zero outside. Both/(r) and F (r) are arbitrary functions 
subject only to the conditions /(0) = and F (a) = 0. 

Putting aside for the present the question of the stability of these simple motions 
or of their resultant, it is clear that if we superpose the two we get another state 
of motion in which we have vortex-filaments in the shape of helices lying on 
concentric cylindric surfaces. The problem to be considered is whether it is possible 
to conceive a similar superposition of two motions in the case of any vortex aggregate 
whose motions are symmetric about an axis. 

There are an infinite number of either ring-shaped vortices, or singly connected 
aggregates (of which Hill's vortex may serve as a type), differing from one another 
in the law of vorticity of the different parts — -the most important being those in 
which the vorticity is uniform. The motions in all these are known in terms of 
the stream function \fj. The value of \jj is however at present only actually known for 
an infinitely thin ring-filament or for a spherical aggregate. 

2. We are to consider two superposed motions. The one component is in meridian 
planes through an axis and can be defined in terms of the stream -function $* The 

■* Throughout ty is taken as the total flow up through the circle whose radius is p. In other words the 

velocity perpendicular to ds is -- ¥ , 

Ji7rp da 
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other component is everywhere perpendicular to these meridian planes. The vortex 
aggregates will be moving with rectilinear translation through the fluid with a 
velocity calculable, when the distribution of vortex motion is known, by Helmholtz's 
method. Bring the aggregate to rest by impressing everywhere a velocity equal and 
opposite to the velocity of translation. The motion then consists of a flow up through 
the centre in the direction of previous translation, the fluid then streaming (in this 
most general case) in spirals round a certain circle. The circle may conveniently be 
called the equatorial axis of the aggregate. The line of symmetry through the 
centre in the direction of translation may then be termed the polar axis. Whether 
we deal with ring-shaped or singly connected aggregates, the surfaces xjj will always 
be ring-shaped inside. In fact they are so also at the boundary, for the surface value 
of xjj really consists in the latter case of the outer boundary together with the 
polar axis. 

3. Conceive now the aggregate divided up into a large number of ring-surfaces 
given by values of a parameter xjj differing by dxp, and confine attention to what is 
going on between the two surfaces \p and xjj + d\p. We shall suppose xjj to increase as 
we pass from the outside inwards. Let dn denote the distance at a point between 
the surfaces xjj and xp + dip, dn to be measured also inwards. In the shell considered 
the lines of flow will be spiral, and the vortex-filaments also spirals, as indicated in 
the figure, the thin line P/ representing a line of flow., the thick Pv a vortex-filament, 



Fig. 1. 




and the line Pm a meridian section. Denote the velocity at P by v and the angle it 
makes with the meridian by <f>. Also let o> denote the molecular rotation at P, and ^ 
the angle the filament makes with the meridian — estimated positive when on the 
opposite side of the meridian to v. 

Consider the flow between the two surfaces xjj and xjj + d\p across the " parallel of 
latitude " through P. The total flow must be the same for every parallel. The area 
through which the flow takes' place is Itrpdn, where p is the distance of P from the 



38 PROFESSOR W. . M. HICKS ON VORTEX MOTION. 

polar axis. Hence 2irpv cos (f)dn is. constant over the surface \jt. It must therefore 
be of the form f(\p) difj. So far xjj is only defined as the parameter which determines 
the particular surface. Choose the parameter so that f($) = 1. ^ is then analogous 
to the stream-function in the simple case. It acts in fact as the stream-function for 
the component of velocity v cos <£. Similar reasoning leads to the conclusion that 
cop cos xfln is also of the form /(*/>) d\p, sayjf, dip. Hence 

2-rrpv cos <j>dn = (ii/> . , . ( 1 ) 

27rpco cos )(dn = t /J dxjj , . . . . (2). 

We started with the supposition that the stream-lines and vortex-lines must lie on 
the same surfaces i//. In other words, there must be no component rotation perpen- 
dicular to \jj. This may be expressed in other words by the statement that the 
circulation round any circuit drawn wholly on xfj must vanish. Take for this circuit 
any two parallels of latitude. The condition gives that the flow along one must 
equal the flow along the other. In other words, the flow round a parallel of latitude 
must be the same for all parallels on the same surface xjj. Hence 

2irpv. sin <f> = f . (3) 

where f is a function of \jt. 

Equations 1, 2, 3 give conditions which any motion possible between any two given 
surfaces xjj and xjj + dxjj must satisfy. In our case, however, the motions in the 
separate shells must fit together. We may regard the vortex-filaments as due to the 
velocities in two successive shells, or as due to the different velocities on the inner 
and outer surfaces of the same shell— the velocities on the inner surface of one bailiff 
the same as on the outer of the next succeeding shell. If now w x be any component 
of a filament, and dA the area perpendicular to a> u the value of o> x dk is given by half 
the circulation round dA. Apply this to the two components a> cos x along a meridian 
and a) sin x along a parallel of latitude. As a circuit for co cos x take two parallels 
one on xjj and the other on i// + dxjj. The flow along the first is 2wpv sin cf> and along 
the latter 



Hence 



But by (3) : 



Hence 



7 

2irpv sin <f> + 2rr -~- (pv sin </>) dn. 
2o) cos x • 2rrpdn = .— 2tt --- (vp sin tf>) dn, 



2-n-pv sin <f> = jf. 

d f 

AlTpCO COS X = — 'J* • • • (4). 
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Comparing with (2) it follows that 



J 1 rl.y 



dn ® dn 

or 



/• — i d l 

Jl - ~ 2 dty ' 

We may regard then Eq. (2) as replaced by (4), which includes it as the greater does 
the less. 

For the circuit for m sin ^ take a small circuit formed by a small arc ds of a meridian 
PP' on \jj, the normals (dn) at P, P' and the portion of the meridian arc on ty + dty 
cut off by these normals. The flow along the normals dn is zero. Along ds it is 
v cos <f) ds ; along ds it is 

v cos </> ds 4- — (v cos (j) ds) dn. 

The area of the cross-section of co sin x is dn ds. 
Hence 



d 

2w sin x dn ds = — -_--- (v cos <j> ds) dn. 



But by (1), 



therefore 



V cos <£ 



dn 

1 rf^r . 

27T/3 $W 



47Tftj sin y efe = — — — as 

A dn x p dn 

Since dty/ds "= i// will give any component of velocity m £Ae meridian plane in 
the same way as the ordinary stream-function. 

4, It will often be found advantageous to express ty in terms of curvilinear 
co-ordinates. Denote these by u,'v. Displacements perpendicular to the u will be 
denoted by dn, and to v by dn, to be estimated positive in the directions in which 
u, v respectively increase. 

The differential equation satisfied by ty is found by expressing the circulation round 
a small area bounded by the curves u,u-\- du, v, v + dv. Let co y ( = co sin y) denote 
the rotation at a point of the area. We shall regard this as positive when it goes 
clockwise. The circulation is then 2aj x X area = 2co 1 dn, dn. 

The velocities along PQ, PP' (see fig. 2) are respectively 

1 dty 1 city 

2jrp dn 2rrp dn 

The llows along them are therefore (clockwise) 

1 city 1 city 

— — L- /-/ nn our i __•-. 



1 city 1 city j 

—-- — - dn and + ~r^~ ~r, dn. 
lirp dn 2wp dn 
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Fig. 2* 




O 



Hence the total How round PQQT' is 



d 



( 1 d4r , A , d f 1 dijr \ 

a% \2>7rp dn j dv \2irp an i 



or 



d / 1 dyfr du 



du \2irp du dn 



dn r \ 1 1 d I 1 d^t dv dn\ 7 7 

. — } dudV — ~r~{ 7 , * 77. 7" ^ «V. 

a?;/ dv ylirp dv dn du J 



But this is 2ft)! c£w dn 1 . Hence 



d I 1 di|r eib flfo'* 






d I \ d^jr dv dn 



die \ p du* dn* dv J dv \ p dv dn'* du 



47TCOJ 



dn dn f 



du dv 



4™ sin x 



dn dn' 
du dv 



(5). 



In many cases p + zl = f {u + m), giving dujdn = dv/dn', and the equation 
simplifies to 

i. (i. ^\ + A / L ^\ = _ 47rai AM 2 . 

du \ p du ) dv \ p dv J \du / 

The following cases will be required :— 
(1) Cylindrical co-ordinates, (p, z) $ 



du = dp = c£n 



and 



or 



d 1 \ d%fr\ 
dp \ p dp j 


j , _1_^[ 
' p d# 2 


d?ty 1 
dp 2 p 


rip tfe 2 



dfo; = dz = dn\ 



4:7ro) sin X) 



= — 4^r-p6i sin ^ 



« » # * » 



(6). 
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(2) Polar Co-ordinates. (r,0), 



p = r sin 9 



du = dr = cfe 



cfo = d!$ 



dn' = r d#. 



and 



d 



6a 



r dty\ d 

wmAnwm 



dr 



1 c£\p 
^ \rp d6 



irrro) sin ^ 5 



or 



<^ 2 -^ 1 <^ 2 ^r COt 6 dty 



di 



de* 



r 2 <i# 



4rrpw sin ^ . . 



(3) Spheroids. 



(a) Prolate* Here p -\~ %i ^z X sinh (it + m) ? 



whence 



p = X sinh it cos 1?, 



2; zz: X cosh w- sin 1;. 



The surfaces u 9 v are respectively the ellipses and hyperbolas 



O 9 

_ I _*_ 

sinh 2 ^ - cosh% 



X 2 and 



A 2 



sin 2 ^ cos% 



A , 



O t - * « 



(n 



u increases from at the origin to 00 at an infinite distance ; v increases from — \tr at 
points on the negative part of the axis of z, through for points on the equatorial 
plane to \tt at points on the positive part of the axis of & - 
Again 



dn\ 2 



du 



dp^ 2 

du 



/ /7 <-v \ £ 



B i B fcj illl— 



d. 



du 



(p + ZL ) ^ (p — *0 



du 



. = X 2 cosh (u + m) cosh (u — - m) 
= X 2 (cosh 2 tt — sin 2 ?;). 

Hence the differential equation is (writing C and S for cosh u 9 sinh u),- 



1 d / 1 <tyA 1 . rf / 1 ^ 

cos v du \ S dfo / S dt? I cos v dv 



47rX 3 o> sin ^ (C 2 — sin 2 ^ ) . . (8), 



.(j8) Oblate. Here p + %<> — X cosh (^ + ^) 



p = X cosh w cos v 

Vfe x 2 



^ 



and the differential equation is 



z = X sinh 'U sin v, 



X 2 sinh (w + m) sinh (^ — m) 
X 2 (cosh% — cos%) 5 



L. d ( \ d ±\ 1 rf / 1 ^ N 



+ 



cos v du\G dtoj l ^ \cos 1? rfv 



4rrX a et> sin ^ (C 2 — cos 2 ?;) . . (9), 



V vJJu# 0-2v.v^Xl.» *■ " "'A,. 



G 
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(4) Toroidal Functions.— Here ['Phil Trans./ 188 1, Part III., p. 614J 



U + VI = 


-. p -f- a -f #6 
= log - j 




P = 


sinli ?.6 
cosh w — cos v 


du 

dn ~" 


smli u 
P 


whence 














d /C 


— cos v dty\ 1 
~~S~ du) + ~S~ 




((C- 

\ 


- cos ^) v -~ —■ 

7 a?; / 




47ra — w sm ^, 










=r- — 


4w s 

(C — cos?;) 2 


o) sin x 



(10). 

5. Equations 1, 3, 4, 5 or 6 give the conditions for a possible motion. It is open 
to us to choose xjj arbitrarily. In this case the equations give v, a>, x* & The motion 
is instantaneously possible, but in general it •will at once proceed to change the 
configuration— the motion will not be steady. The application of this theory to 
values of xjj which are already known (Hill's vortex for example) leads to interesting 
results, but the absence of steadiness robs the theory of importance. If we impose 
the condition of steady motion, it is no longer open to us to choose xjj at will. Let 
us then impose this condition. The condition that the motion shall be steady 
involves : — 

(1) xjj must be a surface containing both vortex-lines and stream-lines. 
This is already the case. 

(2) vo) sin (<f> + x) dn must be constant over the surface. 

It must therefore be of the form Fdxjj, where F is a function of xjj. Hence 



or 



7 I 

vo) sin (<£ + x) ^ -^ 7 * (H)' 

tt ft 

Expanding this, and substituting from 1, 3, 4, 7, 

J cfy [ dr 2 ^ i* dff 1 r 2 d<9 j P 5 

dr* "*" r* <Z0* * r 8 <Z0 ~~~ P J cty ' ' ' ' l ; ' 

where jf and F are arbitrary functions of t// 9 Choosing these, equation 12 will give 
the type of t/f.* 

We proceed to apply these general theorems to certain special cases of spherical 
aggregates. In order to exemplify the method employed we will take first the case 
in which there is no secondary spin, the type in which Hill's spherical vortex is the 
simplest case. 

# For another proof of this equation, due to one of the referees, see end of present paper. 
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Section ii. — Aggregates ivith no Secondary Spin (/= 0) and ivith 

Uniform Vorticity* 

6. We begin with the spherical aggregate, the simplest type of which is the Hill's 
vortex. The equation for iff is that given by equation 7, in which co is put kp where 
k is uniform and x = 2 n - It is 

d:^ i d^ cotedf . 

+ — — ■ — — — — = — 47rA;p = — 47rA:r sin" 0, 
«r- ?^ at/ 2 r 1 dO r ' 

in which is measured from the pole to the equator. A particular solution of this 
x s — |-Tr/<:r 4 sin 4 5. 
In 

^-^ , 1 cPyjr cot dyfr 
dr* + 7 J ~d0* r r d0 = ' 

put t// = r w Z w , Z w being a function of only. Then 

d 2 7 r!7 

— •-oottf^ + »(n-l)Z H =0. 



The integral of this is 



Z w = — sin 6 — — 



where P w-1 is a zonal harmonic of degree w — 1. 
Hence the general solution of the equation in \jt is 



* = - i^ 4 «m 4 + 2 ( A n r» + -?- ) Z, 



j,n-i / « 



Since 



1.3.5. . . (2n-l) f »(»-l) 

| cos- (? - ^-L «**-. + . . . 



7?/ I 



the values of Z w are easily found, except for Z x or Z . It is easily found from the 
direct equation in this case that Z, = Z = cos ft The following results are easily 
deduced : — 

Z 3 = sin 2 0, Z 3 = 3 sin 2 6> cos 0, 

Z 4 = f (4 sin 2 - 5 sin 4 0), sin 4 = £ Z 2 - & Z 4 . 

Consider now first the case of a homogeneous spherical aggregate. In this case 
the functions — Z n apply only to the space outside, and Ar n Z n to the space inside. 

Let ^ denote the value of t// inside and \jj 2 outside. Hence 

G 2 
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^ — — Xjrkr* sin" + 2A„r'Z,, 



& 



c 2 



s * — - z 

7 



Let rt denote the radius of the sphere. Along the boundary of the sphere i//, 
and also dxjj^dr = d\jj 2 /dr. Expressing sin 4 # in terms of Z 2 and Z 4? 






*i 



Itt/^Z, + ^tt/^Z, + 2A H r il Z tt . 



The term in - 1 I r ) 7rir ,] Z 4 may be supposed merged in. A 4 v ,4 Z. b and may therefore he 
treated as absent. The conditions 

xjj l = xjl.. 



U'ive 



and for u 



Hence 



and 



v. 



sz: (i. 



Ax == 0, 



'lr 



A,,, = when n > 2, 
B H = when n > 2, 



A 9 a 2 



2 A a 



|-7r/ca 4 



|-7r)fca 8 



I 



3* 



a 



5l 

rr 



>► 



A 2 — ?,rr/m\ I 






4: —.7.^/5 



i// 2 = -/ 5 -7r£ — sin 2 0. 



The velocity along the normal to the aggregate is 



— - ,/> = i 4 5 -^a 2 cow 0. 
Hence the aggregate moves forward through the surrounding fluid with a velocity 

Referred to the aggregate at rest therefore 

xjj l — |-?T^r 2 (a 2 — r 2 ) mrO. 

The cyclic constant (/x) is the circulation taken round a meridian section, up the 
polar axis and down outside. It is the sum of the circulation round the elementary 
areas of which the section is composed. Hence 
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rn 



Thus 



li = % (elementary circulations) = %2ct) dA ■= &22p <iA 



- %2irp clA = - X volume of aggregate 



m& 



7T 



0) 



A^p 



7T/X 

m 



V 



5^ 



which are Hill's results obtained by direct methods. 

7. Heterogeneous Aggregates. — We may, however, superpose on an aggregate such 
as the foregoing other spherical layers of different vorticities. It will be advisable to 
consider first the case where there is one such layer of vorticity determined by 
(say) k\ We may call them dyads. In this outer portion both terms in Ar n and 
~B/r n ~ l can appear. Let ip u i// 2 , \jj denote the stream functions for each part and for 
the surrounding fluid. Then 



^i 



2. 
5 



Trkr 4 Z 2 + SA w r»Z„, 



^ = - | ff *V% + % (K;r n + ~£) Z B , 



* 



B 



y»' 



Let a, b denote the radii of the two spherical surfaces (a > b), and apply the same 
conditions as before to the two surfaces. 

Again all the co-efficients vanish except for n = 2, and there results 



A 2 & 2 



B./ 



f ir/fe&* = A 2 '6 2 + 



2A a 6"— f-VAfe* 



2A,'6 



2 



B./ 



t TT^fr 



f^'6 8 



> 



and 



^ 



"R "R ' 

-^ = A/a 2 + ^~ 

a 2 s 



^ 



■f Tr/c'a 4 



a 



Bo 



at 



ZijCLo Oj 



B/ 



5. 



•I db'a 8 



of 



5 



.y 



The first two give at once 



the last two 



B,' = -A-ir (k - AO & 8 ? 



A 2 ' = fwa 2 !' ; 



also 



B\ 



= Att {(* - k') m .+ /cv.} 

A 2 = ■■-§-* .{Fa 2 -f (F- F) & 2 }. 
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Met + (k - h f ) ¥ 



,2 



o 



jh A I sin 2 0, 



xjj 2 = 2it \ ^k'arr 2 + ^ (A - F) 
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PV 4 !> sin 2 ft 



* 



— _ „ glir 

Id r 



The normal velocity at the outer boundary is 



(when r = a) — ~ x % — — '— cos ft 



2-7775 rdO 



cv 



The outer boundary therefore progresses unchanged with velocity of translation 



v — 15 



(k - V) ¥ ~h Va? 



a" 



Bring the outer boundary to rest by impressing on every part of the fluid a 
velocity equal and opposite to this, i.e., adding to the stream-functions a term 



T~5 ^ 



(k — k f ) Jf + k'a h „ . 9 A 

- — ™-rsm z &. 



«° 



The relative motions are then given by 



I7T 



^ = -j \± ( 5&'a 2 + 5 (A - *') ft 2 - 2&V - 2 (& - A') »■ 



<r 



2tt 

I 



y UV-^+Ki-F) 5 



2P 



»'- 



r 2 sin 2 



l//r 



' I' (a 2 - r 2 ) f* + I (A - &') ~ $: (a 3 - r 8 ) 



<rY 



► sin 2 ft 



&?* 2 1 r 2 sin 2 



If, however, the motion is to be steady, the inner sphere must now be at i^est, 
that is \\ii = when r = b. We get, therefore, the following necessary relation 
between k 9 k\ a, b $ 



Jc'a 2 - Job 2 + i(k- h') 5 



This may be written 



^) b 2 = 0. 
a* J 



2b 2 k (a 8 - V) + k f {Sa s (a 2 - b 2 ) - 2& 2 (a 3 - S 3 )} = 0. 

Both the expressions in the brackets are positive, hence kjk' must be negative or 
the rotations in opposite directions in the two portions. 
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Denote the cyclic constants of the inner and outer portions by [x ly jut 2 - As before, 
we see that they are respectively 



h 



IT 



X vol. 



That is 



//,! 



AiO 



IT 

m'k' 



— — o U IV, 



4 



7T 



| ( a 3 _ & 3) ^ 



Substituting for k 9 k f in terms of ^ 1? /x 2 



6 2 



a" 



& 5 



v^i^ij + ^^.p) 



The result is that a double aggregate is possible. If, however, the size is given 
the ratio of the vorticities must have a special value, and vice versd. In terms 
of the radii it may be shown that 



V = — 



4 (hj-V) b 2 (a - b) {2a? 4-_3& 8 _H- 4o%_+ 6a6 2 } 

45a 3 ~+ 6) " 



Three cases specially invite attention, (1) equal volumes, (2) both parts made of 
similar matter, i.e., vorticities equal, and (3) equal cyclic constants. 



Case i. — Here a 9 = 2b l 

h z 



w 



i 



66» (a 2 - 6 2 ) - 26 5 



3 x 2* 3 - 4 






76220 



|- nearly. 



Case ii. — h' 



k 



3a* (a 2 - b 2 ) - 46 2 (a 8 - & 8 ) = 0, 



Put alb = x, we get 



3x 4 + 3x z — 4.x 2 — 4x — 4: 
This has three negative roots ; the positive one is 



0. 



x = 1*3283 or -^ = f nearly. 



a 
J 



Case iii. — /x = — \xl or 



h 



a s — ¥' 






h - fc' 



&° 



whence 



a *b* _ fe2 ^3 _ 6 s j + i a a 5 _ 



W 



a" 



>\ 



>2 — 



0, 



2a 3 + V — 3a?b = (a - b) (2a 2 — a& - b 2 ) = (a — &) 2 (2a + &) = 0. 



Equal circulations are therefore impossible. 
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8. Poly ads.— "Passing on now to the consideration of any number of layers, let the 
radii of the spherical boundaries from the inside outwards he denoted by a l9 a 2 . . - ci n ; 
the vorticities by h u k 2 . . . k, M and the stream-functions by xjj^ \jj 2 . . . \}i n and \jj. Then 



- ( 



J li ^ XI. i 



A,r 2 — i hr* } sin 2 $ 1 



\ B 1 

y p = 27r« A^r 2 + ~~~- — I /^ 9 r 4 j^sin 2 </ , 



I 



ik +1 = 2ir -^ sin 2 

Applying the conditions of continuity at the pfh boundary , there results 



B 



T> 



with 



Clnj lip 

9 A r* 2 — — — - 1 ^y 4 — - 9, A a 2 — -^ — ^ & a 4 



-m^ -f I 



Adding 
Similarly 



. p+1 — A p = |- (^ +1 — k p ) a 2 p with A n+1 



2 
'jp + 1 '" ~ -^ — " ~ 1 



t (K+i ~~ K) a n wit 'h B i = 0- 



Clearly the A ? s evolve from the outside, the B ? s from inside. 

Write 

\(k p ~k p+1 ) = \ r 

Then 

J\.p *°~™ Apij ^^ A.pCO p Wit II A w ^_j ^^> U ? 

B p+1 — Bp = f X p a p with Bj = 0. 



Hence 

4 — w *v w ^ — r , ■, ** w 



XX -—— • ^rp lXpVVrp ^ -*-^'}3 """""" Ri ^1 "*""f)^% 



Thus the i|/ are completely determined . 

For steadiness of motion it is necessary that the translatory velocity of the different 
boundaries be the same. This is obtained if the velocities of the inner and outer 
boundaries of each layer are equal. 

Hence we get n — 1 equations (p = 2 to n) 

B B 

i V = A,, + -f - i hjfil = A p + -^- - i ^a|_„ 



or 






or 



9 9 






# ^«. 
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If the volumes of all the layers are equal, 



Hence 



or 



Now 



xxeiioe 



at 



p 



al_ x = a\ and a* = pa\. 



B, 



p 



; Kv (P ~ 1) {P m - (P ~ 1) 2/3 } «? 



p(p-i){p^~( P -iy*}i< P 



2 {( P - i) s %_i + ■ • • - + M- 



^>-i — 3 \fop-l ftp)' 



Up -i)p" is - (p + i)(p~if^ k p =~ Up -if%^~- 2{(p-2f'\_ 2 + . . . , +\,\ 
= - 1 (( (p - 1) 5/3 - (p - 2) 5/3 ) h-i + ••.. + (2 6/3 - i)k + *>}, 

or subtracting two consecutive equations 

{p w _ (p + i) (p _ i)^} ^ = - {(p ~ 2 r - ( P - f) (p - if 8 } Vi. 

Thus the & can be determined in order from the inside. The peculiarity is that the 
process can stop at any point. That is that if we have two poly-ads, with m and n 
layers respectively (m > n) then the first n layers in the first will be precisely similar 
to those in the second. The values are 



and when p is large 



ho 






3 x 2 2/3 - 4 
+ 1-4717*! 
— l-5866*i 



w± 



h» 



'£?)-!• 



1-3120*! 



As another example, take the case where the layers are formed of the same 
material, i.e., the vorticities alternately equal and opposite. Then k p = ( — Y" 1 ^ 



^P — 3 I^P — 3 ^1 ( ) UWL A n — g- fC n — - ( y 3 ^i 



3 ,XJ % 

T~5 \ a i "~" a 2 4" • • • ♦ + 

Let x p denote the ratio a p+l /a p . 
These values are then given by 



P-2/Y. 5 



^jp-l) 5 \ / ^2^-1 



2 9 

CC p ttp~i 



Ju 



9 

3 P 



P 



X 



p 



l 



|U 



«£" 






1 



2)-l 



/>•». 
(//. 



tyvi 



.... 



and may be found in succession. The equations are, if b p denote 1 



x p v^p \ 1) — 3 ®P \ X P i x p ~r !)• 



vol. oxen. — A. 



H 



a? 



I • » a 



'# — 1 
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In which it is clear that 



If b p = -J, the equation is 



6=1- J *- 1 



X % (X + 1) = f (# 2 + # + 1): 

the positive root of which is a? = .1. In this case 



frp+i 



1 — --2- 

1 



1 

2 



If ever ^ is nearly -J = 1 -f a (say), a? p is nearly I = 1 + £ (say). Then, regarding 
a and f of same order 



( I + 3f + 3f ) (2 + f) = f (3 + 8f + ^) (I + 2«) 



£ = £* + l*f--£f 



4 



7T a " ™ 5" a 



(i-A«) 



and 



0«+i "— • I """" 



i + « 



{H{« 



75 M i 



l-(±+*)(l-4* + W* 2 ) = * + «- 



15 a# 



Hence 6 P continually converges to J and the value of x p to 1 asp increases. 
The first seven values are. 

x. 2 = 1-1840, 
a 8 = 1*1284, 

a;4=l-0987, 
05,= 1-0802, 
jc 8 = 1-0674, 
057=1-0580, 



6,,= 


•7582. 


6,,= 


•6741. 


6 4 = 


'6315. 


6 6 = 


•6056. 


&„ = 


•5882. 


6 r = 


•5753. 


& 8 = 


■5660. 



The succeeding values will be given to four figures by the foregoing approximations. 
The velocities of translation of the series of aggregates are 



Monad 


v,= 


1 5 ^l a l 


• 


■ — 


V,. 


Dyad 


Vo = 


4(7- 


5a-?) V, 


— r: 


— -9110 V,. 


Triad 


v,= 


4(7- 


lOa-f + 5ajfa^) Vj 


- — 


+ -8615 Vl 


4 -ad 


v 4 = 


4(7- 


lOxj + lOxfxi — 5x{xlxi) Vj 


—- 


— -8282 V,. 


5 -ad 


v 6 = 


4(7- 


I Ut^j "~|" IUt^jtA/2 "*""* I \)0ui9&2&3 "j odu^Xj^Jb^Ju^\ \ i 


:zz 


•8023 V,. 


&c. 


v 6 
v 7 






= 


- -7833 Vj. 
•7618 Vj. 




v 8 






^= 


- -7462 V,. 
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9. The form of the stream lines for a monad aggregate have been delineated by 
Hill. The general form of the stream lines for a poly-ad is obvious, and there is no 
special reason for drawing them accurately at present. It will be well, however, 
to determine the position of the equatorial axes, for the particular case of homo- 
geneous poly-ads, that is in which h p = { — ) p h x . 

The condition at an equatorial axis is that 

7~~~ -~~r = when 6 = — , 

2irp ar I 

in which i// denotes the stream-function referred to the boundary at rest. Applying 
this to the p-th layer in an n-ad 



\jj p = 2rr 



B B r* 



' ^ •» 



The equation for the equatorial axis is therefore 

- B, (-~ + |r) - 1 V' 3 + 1 V<2 = o, 

or 

r 5 (a» - ^_i) - I (a* - ajL,) '/' 3 - | ^jLi (4 - <-0 = 0. 

This may be written 



Now 



therefore 



and 



n£ __ JL f£=l_* ± r ,2 ,,3 _ 1 ggrl J" -.3 „5 __ n 

^ 2> _i (a^>-i — 1) = s : ^1 (^)-i — 1), 

T r> — ~ 1 



For a monad 



2 
2 Cfc o& 



Q y zzzz o* — •— 



for a dyad 



v' 



2 



&i = 1 r 5 — | afr 3 — I w{ = r = 1-1720^. 
Beyond dyads 

r = Op^x nearly = (1 + £) a 2)-1 . 

JL lit? IX 

<i { 1 ~ i (I &,-i + 1) ~ s L &„-i} 4- <,f {5 - | (| 6„ ; + 1)} - 0. 
Now b p „ 1 is nearly |- 
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The distance from the inner layer is therefore 



From the outer it is 



But (p. 50) 



therefore 



2 



Tp <X>p-l — sjp-l^'p- 



Tv—i"-'v—l' 



\Alrp "*"*"* 7 



1 J 



Clp-.\ 



5Jp-l a p~i 



Ratio = 



2. f 

5 J V-l 



Jj 3?-1 X 5 J3?-l 



X„i -1 + sA 



'2)-l 



1 5Jp-l9 



Ratio 



or the equatorial axis, with increasing number of layers, tends to bisect the distance 
between the two boundaries of the layer. 

10. Energy.— The energy within any region is 



E=| 



fJ-J^ + ^J}^*, 



J(27rp) 2 l\^7 



the integral extending within the boundary of the region. By the ordinary method 
this is reduced to the form 



— — '■— ds + o)\Ij dp dz. 
:tt J p an J J 



Since 



4z7T J /5 

*L l L <*t\ i A / 1 ^ 

dp \ p dp J dz \ p dz 



Arret), 



If the boundary be infinite and the fluid at rest then the first integral is zero, 
and 



E^ 



O) 



xjj dp dz. 



The integral extending only to spaces which contain rotational motion, 
motion is of uniform vorticity co — hp, and 



If the 



-R = k 



j p^dp 



CbZe 



In the cases here considered ip is of the form f(r) sin 2 0, and 



E = 2k 



fr*f(r) 

JO 



dr sin 8 6 dd 



4 L 
3 K 



r 2 f(r) dr. 



In the case of a poly-ad f(r) is different for the various layers, and 



E = f 



h 



% 







*a 



r % fi (r) dr + k 2 \ r 2 / 2 (r) dr + . . 



% 
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We work out the case for a dyad aggregate, in which k 2 = — k u 



5 ^V } 



A(r) = 2rr{(§a\-$a%)k 1 r*- 

A (r) = 2ir { - i k ia y + -A- half 7- + | ^r 4 } , 



and 



E 



8w*f 



o 
O 



(ts- (2a] — of) a? - ^ a? + fV «1 (a| 



a?) 



1 5" ^l (^2 ^1/ ~~ 5" • 7 1^2 """" Ml) $ 



15 



{ 3 a i 3 a 2 a l ~f* "3TT a 2 j 



32^ 
45 



(ai — a^af + | a?) 



1 OK OuTTK^ „ 

= l*y4o X , - ~ ce x . 

If the two parts had been single monads their combined energy (when far apart) 
would have been 



E = 1-534 X 



'/.2 



45 x 7 



a{, 



The energy when combined is therefore greater than when they are separate. 

11. It may not be out of place to make a short digression here as to the relation 
of a Hill's vortex to the vortex rings which have been investigated in previous parts 
of these researches. As is known the translation velocity of an ordinary ring 
decreases as the energy increases, and formulae are given in a former paper* whereby 
those quantities can be calculated for comparatively thick rings up to R/r = 4 with 
considerable accuracy, and possibly further. Here R is the radius of the equatorial 
axis and r the mean radius of the section of the ring. Refer all measurements to the 
spherical form, and let c denote its radius, V its velocity of translation, and E its 
energy. Take now a ring of the same volume and circulation as the sphere* and let 
V" and E denote its translation velocity and energy. We get the following value of 
E/E 0j V/Vo for different apertures. 



R 


R 


r 


V 


E 


r 


c 


• 

c 


V 4 


E 


100 


• » 


• • 


•199 


176 


50 


8-09 


•162 


•282 


95 


10 


2-77 


•277 


•593 


20-8 


5 


1-745 


•349 


•784 


10-25 


4 


1-500 ; 


•375 


•856 


8 


3 


1-239 


-413 


•946 


6 



* "Researches in tlie Theory of Vortex Rings," Part II., p. 757, ' Phil. Trans.,' 1885, Part II. 
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These numbers are graphically represented in fig, 1, Plate 1, where the abscissae 
give E/E and the ordinates V/V . Dotted lines refer to points where calculation 
cannot be applied. On the same figure are placed outlines of the aggregates drawn 
to scale. Two things at once strike the eye. First, that the spherical aggregate 
evidently lies on the E .V curve of the rings, belongs, in fact, to the same family ; 
and, secondly, that the variation of V with the energy is small over a very large 
range. The shape and nature of the aggregate when the energy is nearly that of the 
spherical form have not yet been determined. It is probable that as the energy 
diminishes the form lengthens along the polar axis, until when the energy is very 
small it becomes a long, thin, cylindrical aggregate. When this is so long that the 
end portions form only a small portion of the whole, it is possible to obtain an 
approximation to the energy, for when very long the fluid outside will be very 
nearly at rest (as in case of force outside a long helix). The velocity of propagation 
will then be the velocity at the axis. Let a be the radius of the cylinder, I its 
length. Then 

]ri 2 _ si 



Id = 4c 3 , 



Again, if V denote the velocity along the axis, the velocity outside is zero, and the 
variation at the ends only a small part of the whole. Hence the circulation is 
given by 

Again let v denote the velocity at a distance r from the axis. Take a small 
rectangular circuit, b parallel to the axis, one inside distant r from the axis, the other 
outside. The circulation round this is bv. But it is also the value 2,codA. taken 
over the area of the rectangle. 

Therefore 

k k 

bv = Jc£2r dA = — (volume) = . brr (a 2 — r 2 ), 

rrr ^ TT 

v = k (d — r 2 ), fji = Ikd 2 ; 
therefore 

V = — 1 :; 

I \ a 1 i 



Energy in E ~ 2nr . I dr . -^v 2 

Jo 



*> \ o 



¥ r. i 1 - $) d c") 



TTfjira 2, [j? 9 n c 



61 ' U l 



. m = fTTfJL 



x w y2 _ | ^3 Y 2 t 
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It is thus the same as a mass of one-third its own mass moving with its velocity of 
translation. Now 

E = 



therefore 



E 



therefore 



^ 



h V c > 


V — i 6 - 

5c 


B5 _ 6 _ 


V he 


9 /2 ' 

V 


v. / ' 


E _ 7 


/ V \ 2 

1 • 



K ~ 45 \ V, 



i^ 



This only holds, however, when V/V is small. It is a small part of a parabola in 
the figure touching the axis of E/E . 

12. Spheroidal Aggregates.— As is known from Hill's investigations, the spheroid, 
although an instantaneously possible form, is not steady. It proceeds at once to 
change its shape into a non-spheroidal one. It seems, however, advisable to give the 
general outline of the method as adopted in this paper and as applied to the 
spheroids, in order to investigate whether by superposing a second or third layer it 
may be possible to obtain a steady form. 

The functions involved and the differential equation for xjj are given in Eqs. 8, 9. 
Writing C for coshtt and S for sinh u> the differential equation in \p is 

— y (± ?) + I f (— Tl = - 4^-ibX* S cos v (C 2 - sin 2 v), 

cos v du \ S du j b dv \cos v dv j K ' 

since 

o) = hp = /cXS cos v. 

As in the former case, a particular integral is 

xjj -=z — — - kp A = — 1 7T^X 4 S 4 cos 4 V. 

It remains to integrate 

1 d 1 1 dyjA J_ d^ I 1 ^ 



+ — - ; -—= 0. 

cos v du\ $ du) S efo \eos t> cfo / 

This can be satisfied by writing \\f = 2X m Z OT where X and Z are functions respec- 
tively of u and v only, and 

d I \ dZ\ mZ 



d-v \cos -v dv) cos ?; 

^ /J_^X\ __ mX 
<#w \ S du I ~ S 



)►. 



w- being any constant. These equations are 
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CvZt cLii r7 

-r~ + tan v — - + m/j = 

dv av 



-\ 



d 2 X J 1 dX 

v^ ~ - cotn u — 
du u du 



> 



mX = 



As will be seen later m must be of the form n(n-— I), n being any integer. 



7T 

Writing for a moment v = -r- 



#, the equation in Z becomes 



whence 



Therefore 



cWL 
dff 1 



dZ 
d0 



cot 6^ + n(n — 1) Z = 0, 



Z 



w 



/7P 
Sin V — r-r" 

dd 



Z 9 = sin 2 # = cos 2 v 



Z, 



6 sin 2 5 
6 cos 2 ^ 



- 1 / sin 4 6 
-£• cos 4 i; 



cos v = 



4. 7 _2_ 7" 

5^2"" 15 ^4° 



To determine X, we proceed by the same analogy to put 






dP 

du 



Then 

(FX 
du 2 



dX 



coth u— ~~ n(n — 1 ) X = S 



die 



du 



d^ , t . dP 
-7-7 + cotn u — 

air <m 



n{n — 1) P 



?. 



If then P denote a zonal harmonic with imaginary argument and of order n — 1, 
the right hand of the above vanishes, and the value of X is a solution. That is 



dV 

du 



Now we have 



P 



1.3...(2ra~ 1) 



w 



" s 



%I 



c „ 



n (n — 1) 
2 (2fl,~^T) 






Hence 



Q2 

2 - — O . 



"V — - £?Q2 1 1 5 Q4 04 __ 2 *Y 



5 ^2 



This set of solutions gives values finite and continuous at all points inside a given 
ellipse of the family, but infinitely large at an infinite distance. Let Y denote the 
second, integral of the equation. Then, in the usual way, it may be shown that 

rs 

J X 
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whence it is easy to prove that 






V — Q 2 — — 1 ft 2 1r>rr Z- -_ 1 P 



C 4- 1 

J. 4 = "24 -A_4 log p ~ — -gij: Lv ( LoO — - loj. 



ith these values the particular integral is 



2 
2 



5 - tt&X 4 (2X 2 - I X 4 ) (2Z 2 - i Z 4 ). 



1 The terms in X 2 Z 2 , X 4 Z 4 may be supposed merged in the general solution. We 
mav then write 

^ = ( A 2 X 2 — - 7 ~ 5 - 7riX 4 X 4 ) Z 2 -f- ( A 4 X 4 — Y5 7r^X 4 X 2 ) Z 4? 

^2 = B 2 Y 2 Z 2 + B 4 Y 4 Z 4 . 

From these it is easy to deduce the values of A 2 , etc., for a single free aggregate 2 
by applying the conditions \fj 1 = \fj 2 and dxpjdu = d\jj 2 /du at the surface. It is 
unnecessary to do this, as from Hill's work we know that it is not steady. 

The case of motion inside a rigid spheroidal boundary is also given by Hill.* 
The solution follows immediately by impressing the condition ^ = when u = u. 

Hence 



A 2 — ~f~5 Trfch — 

X 2 



A, = -A- irkX 



4 _7*\4A2 

-A. 4 



7 5 



? 



where thick type denotes values at the surface, and 

Y Y Y Y YY Yy 

yi — 75 7r/1/A " xf £**£ 75 7r/(/A • ~~ y ^4? 

A 2 A 4 



which easily reduces to 



^ = J?^L (S 2 - S 2 ) S 2 (S 2 + cos 2 v) cosV 
4 -f- 5S" 



The total circulation is f- &\ 3 CS 2 . 
The equatorial axis is given by 



._= 0, when v = 0. 
du 



That is by the equation 



2S 2 S- 4S 3 = 0, orS~ 4|S. 



* < Phil. Trans.' Part II., 1884, p. 403. 
VOL. CXCIL— A. I 
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The equatorial axis therefore lies in the equatorial section in a similar position to 
that for a sphere. 

13. Dyad Spheroids,— Poly-ad spheroids clearly occur in the same way as for 
spheres ; they are, however, also unsteady. It will be sufficient merely to indicate 
the steps of the proof. 

Let u = v! and u = u" denote the two boundaries. \p will involve terms in Z 2 and 
Z 4 . By applying the surface conditions in the same way as for the spheres to both 
sets of terms independently, the coefficients are determined, whilst the condition that 
the internal interface has the same translational velocity as the outer gives for Z 4 an 
equation which xi and u" must satisfy. This is 



w T h 



ere 



C'S' 1 ' 
S" 


M" 


— 


C"S" J 
S' 


M' 


+ 


(TS /;i - 
6 - V 6 


C'S" 
S' 2 






M 




du 




• Y 4 


du 


? 



o, 



and the dashed letters refer to values at the outer and inner boundaries v!\ u. 
The same applied to the Z 2 terms give 

~S" N" - ^r N' + (C"S" 4 - C'S' 4 ) ^ = 0, 



where 






The existence of steadily-moving spheroids depends on the possibility of finding 
values o£u\ n f to satisfy these two equations. 
It is easy to show that 

6M + N = |-S(25S 2 + 14). 

Hence, adding 6 times the first equation to the second, there results an equation free 
of logarithmic terms and which can easily be reduced to 

S' 3 (S" 8 - S /2 ) _ 2C 

C"S /A1 ~ C'S /4 = 50' 2 "'~i ' 

Putting C = y, G" = x, the factor (x — yf divides out, and the equation may 
be put in the form 

2y (x s + 2x 2 y + 3xy 2 — 2x) + (y 2 — 1) (3^/ 2 + 1) = 0. 

Now x > y > 1. Hence 3^t/ 2 — 2x^xy 2 + 2$ (y 2 — 1) is positive. The expres- 
sion on the left is therefore always positive and no suitable values of x 9 y satisfy 
the equation. A prolate spheroidal dyad is therefore not steady. 

The condition for the oblate spheroid can be found by writing S-y/ ~~™ 1 for C. It 
can be shown that this also has no suitable root, 
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Section iu.—Gyrostatic Aggregates. 

14. Passing on now to the consideration of the more general problem where a 
secondary spin exists, the simplest case is that in which in equation (12) both P and 
fclfjclxfj are uniform. 

Suppose 

/J£ = A, or /=V(2A*). 
The differential equation in \b is now 

d 2 y]r . 1 cP"\lr cot d^fr 

a particular integral of which is 

and the general integral is the same as that considered in the pi^evious section, viz. : 

T / 

It will however not be found possible to satisfy the boundary conditions unless the 
term \jj = A x r be introduced. This term, as well as that in |Ar 2 } makes the motion 
discontinuous at the polar axis. However, we will suppose for the moment this 
portion of space excluded, and see later if it is possible to do so. The stream- 
functions are then,- — inside 

^ _ _ „y F _ i A ? ' 2 + A l r + t 2 A n r'%, 

outside 

B 4 



'M-5 



T2 * — ^ n-l ^ 



and p 4 can be replaced as before by f r 4 Z 2 . 

Applying the conditions ^ ~~ i// 2 and dxftjdr = d\\f 2 jch% when r = a it is easy to 
deduce that 

i/f, = - £A (a - r) 2 - |-7r 2 Fr 4 Z 2 + f 7r 2 aWZ 2 , 



fc = - X V 2 F ^ Z 2 . 



The velocity normal to the sphere is 

1 d^f 



-j^g- TrFa 2 cos 0. 



r=a 



27T/) n£0_ 

That is, the sphere progresses bodily with a velocity given by 

V = -& 7ra 2 F. 

I 2 
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Impress — V on every point, that is, deduct yj 7r 2 a 2 Fr 2 Z 2 . Then the stream- 
function referred to the boundary is 

i//j = -|A(a- rf + 1 7r 2 Fr 2 (a 2 — r 2 ) Z 2 . 

At the outer boundary \jj = 0. If we trace the stream-line t// = 0, it is seen that 
it consists of the circle r = a and the curve 

•|A (a.— r) = |?r 2 Fr 2 (a + r) sin 2 ft 

This passes through the poles (r = a, 6 = 0) and touches the circle there, Hence 
the space between this and the outer boundary does not contain, the polar axis. The 
motion given by \ft is therefore finite and continuous there, The space inside it 
must be excluded as giving a motion not possible— or rather, a motion due to sources 
and sinks on the polar axis. We shall suppose it excluded by replacing the fluid by 
a solid nucleus of the shape required. 

The radius of an equatorial axis is given by dxjj/dr = when 9 = tt/2, or by 



A (a - r) + %ir 2 Fr (a 2 - 2r z ) = 0, 

5A 
In this write via = x and —-77—* = b. Then 

x A + (b — |) x — & = . , (13). 

This has one root between and 1. The other roots must either be both imaginary , 
or, if real, one at least must be negative, since the coefficient of x? is zero. As, 
further, x = — 00 and x = both make the expression on the left of the same sign, 
both these roots must be negative. Hence there is one and only one root between 
and I. That is, there is only one equatorial axis. 

In the special case b — \> the radius of the equatorial axis is a . 2~* ~ '7937a. 
For this curve 



r l 



%b v ■= -i-A < --j (a 2 — r 2 ) Z 2 — (« — r) 2 

The curves are drawn in fig. 1, Plate 2, for values of 2t/f/Aa 2 = — 9 1, 0, + '1- 
The value at the equatorial axis is '397. The value (— °1) is drawn to show how the 
discontinuity enters. 

The velocity along a parallel of latitude is given by the equation 

2irpv sin <fi = / — v /(2Ai^). 

This is zero at the surface and on the spindle-shaped nucleus, and increases to a 
maximum at the equatorial axis. The secondary cyclic constant is the circulation 
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round the two circles (1) the equator of the sphere, and (2) the equatorial axis. It 
is therefore given by 

or 

where x is the root of equation (13). 

On account of the artificial nature of the internal nucleus the farther discussion 
of this case is scarcely called for. We pass on, therefore, to the more important 

case — the next simplest one — -in which F is uniform, but the second, terms varies 
as \fi. 

1 5. Case fy-cc xfj, — Here also f varies as i/r. 

Write f = —\Jj where a is a length, which may be taken to be the radius of the 

8ttW 2 
sphere, and X is a pure number. Also write F = -==— where V is a velocity. Then 

the equation in \jj is 

cP^ f 1 d*f cot 0dyfr _ /> 2 Tr X 2 

a 



dr 2 * r % d6 2 r % dd ~ ~ a 2 ™ % ^* 



A particular integral is — — p 2 and the general integral depends on 



V 



!** + I s d6 % V d0 + ~^ Y ~ 



In this put \fj = J n Z n where Z n is the function of 6 already discussed (§ 6) and J n is 
a function of r only. Then 

<^ 2 J» \n{n —1) X 21 



dr % r 2 a 



J, = 0. 



Ji/\/r is therefore a Bessel's function of order n — |-, which can, as is known, be 
expressed in finite form involving circular functions. In what immediately follows, 
the values of J 2 will alone be required. The equation is, writing x for r/a, and 
dropping the subscript 2, 



Ci"u I Li \ 9 \ T 

"TT -— ( 7^ ~ X , i J 



0. 



If J and Y denote the two integrals 



k# 



Jbi.ll. /V<a/ . 

= — COS AX, 

Xx 

= - -; - + sin Xx, 

Xx 



uZ 
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or, in more general terms, 



Integral = C 



sin {a -f \x) 



<, - 



\x 



cos (a + Xa?) 



where C and a are arbitrary constants. 

J and Y may be expressed in infinite convergent series. Thus 



-j- , x SID 'W 

J (y) =E — -- — cos i/ 



Y(y) 



also, 



and 



3 y 



1 9 



-ij . »J . D 



/+...(-) 



ft + 1 



2n 



(2n + 1) ! 



y 



:n 



10 



+ ...+(-) 



m 






2/ 



2m 



2m + 3 (2m + 1) ! 



4- 

I * « « 



cos 2/ 

y 



+ sin i/ 



1 



y 



+ iy + • • • + ( — ) n+1 



272, — 1 

(2w)! 



r 1 "" 1 + . 



2n — 1 ' 



dY (y) 

cly 



sin y 



cos 2/ 



Y 



AT 



•J 



2/ 



*> 



-*r-'y* 



sin. ?/ 

.9 



r 



cos y 

"7" 



>■ 



i 



* # « 



(14) 



(15), 



(16). 



Clearly the functions J refer only to space excluding infinity ; Y to space excluding 
the origin. 

16. For the problem in question the stream-functions are, therefore, 



inside, 
outside, 



*i = 



- 7 r 2 sin 2 6 + %AJ n Z n , 
X~ 



13, 



1 f. — ^ _±?2L £ 



Applying the surface conditions that when x = 1, ^ = i/% and c^/db = d\jj 2 /dx, it 
follows that when 



n > 2, A w = B w = 0. 



w 



hen 



?i 



V 



R 



X 2 


ct 


M'MI liw— /" * 


■2° — 


a 


J 


2V 
X 2 


a 2 


-LA 


dS' 

dm 


•*t/*&K&8B> 


J.Jft 



where J' and d¥jdx mean the values of J and dJ/dx when x = 1, that Is 
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63 



T/ sin X x 

J = — cos A, 

X 



CltJ \ • \ Sill A . x \ • \ T/ 

— ■ = X Sill X ~- + COS A = A Sill A — • J , 

clx fa 



the two equations for A 2? B 2 give 



XX'; 



3V 



X s sin X 



a 2 , 



B, 



V« 8 / 3 J' 



X a VXsinA 



The aggregate moves through the fluid with a velocity of translation given by 



u 



2B, 



V / 3J' 



27ra s 7rX 2 \\ sin \ 



- 1 . 



By its formation the above value of xfj satisfies all the equations of condition except 

that in those equations \fj is the velocity- function referred to fixed axes. Here it is 

not — it represents the motion referred to the instantaneous position of the sphere. 

It is, therefore, not directly applicable unless the velocity of translation given by it 

vanishes, that is, unless 

J' — ^X sin X = 0. 

If X be a root of this equation we get a steady motion of a vortex aggregate, at 
rest in the surrounding fluid. 

If we, however, take the above general function, it gives a velocity of translation 



U 



Y 



3J' 



7T A- 2 \ X sili X 



• • * • l JL f I. 



Bring the aggregate to rest by impressing a velocity™ U on the whole fluid — that 
is, add to the stream-function a term — - 7rUp 2 = — 7rUa 2 x 2 sin 2 0. 
We get a new value of \f), referred to axes remaining fixed, viz., 



* 



3W 2 
X 8 sin X 



(J - x'J') sin 2 9. 



X 



Take this value of \jj, and put f = -\fi. Then equations (I, 3, 4, 7) become 

, 1 clylr 

vp cos d> = ~ - 

r ^ 2?r dn 



vp sin <f) 



cop cos x 



o)p sm x 



9 



ira 



X dyjr 



4c r ira dn 



3V 



o 



4tt\ sin X 



J sin 2 0. 
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These give v, w 7 <£, % 

Now substitute in vco sin <£ + ^ The result is that 



vo) sin (6 + x) dn = 



3V 

8tt% 2 'X sin X 



J'.tty 



so that the motion given by the new xjj is a steady one, There exist 3 therefore, 
systems travelling through the fluid with velocities given by (17) and with a steady 
motion. The system given by J' = i|\ sin X is contained as a special case, 

17. There are two circulations to be considered. That along a circuit up the polar 
axis and down over the surface of the sphere, and that due to the motion round the 
polar axis* Call them respectively the primary and secondary cyclic constants, and 
denote them by /x, v, 

1 Hr } j M i o P /2 f JL. !b£ 

%irp dr 



p 



=,f 



o [ 2ttp r d6 



\ dr + 2 ["'{ - ~ 

J 0=0 Jo L ^' 



add. 



r~a 



In finding this the term x 2 3 f sin 2 may be omitted as giving no circulation, and 
we may take 



where 



Now 



therefore, 



*= 



3W 



X 8 sin X 

3Va 2 



J sin 2 6 



ix = 



7rX s sin X 



2 h dr 



(W f "sin cl6 

Jo 



3Va 



...... j 



tfX 2 sin X 



Jo 



21 — dy 



r 



dr 
dT 1 



y = 



Xr 



a 



\j<*y = 



2/ J ^ <% 



* = - 7 + i(-f - y) **• 



2 f.7 * = ~ L7 L+L ir* 



Also, J (y) is of the order y 2 when y is small, therefore, 



and 



2 iy d y=~i + 6A - 



SVa 



M- = 



?rX 2 sin X 



(S^X — sin X). 



If we replace V as a constant of the motion by /*,, 



*== 



<?r^# 



X (SiX — sin X) 



(J - a> 2 J') sin 2 $ 8 
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Before discussing the value of v it will be well to get some general idea of the 
nature of the motions. One of the most striking peculiarities of these aggregates 
is the quasi -periodicity of type as X increases from to infinity. The best way to 
illustrate this is to use a graphical construction. Now 

xjjaz (J(X^) -a?J(\)}. 

In fig. 2, Plate 1, the curve y = J (X) is drawn. P 1 corresponds to a given type 
(X) of aggregate. A parabola is drawn with vertex at O and passing through P 1# 
Represent any abscissa to the left of X (or of P x ) by Xx, where x<l. Then the 
differences of ordinates between the curve and the parabola up to P represent 

J (Xx) - x 2 J (X). 

It is clear from the figure that, in the position P 1? this function never vanishes for 
x<l. In the second position, P 2 , however, the parabola intersects the curve at 
another pointy. For this point (suppose x = x ) \\s vanishes for all values of 0, and 
the corresponding current sheet is a sphere internal to the boundary. The aggregate 
consists of two portions with independent motions. The primary circulations are in 
opposite directions, and there will be two equatorial axes. So, as P moves on 
along the curve, i.e., as X increases, we get families of aggregates with three, four, &c, 
layers, and a corresponding number of equatorial axes. We shall denote any transition 
value of X by X 2 . Each layer will have its own secondary circulation, given by the 
circulation round the double circuit formed by its equatorial axis, and an equator on 
its boundary. 

Now the secondary spin velocity is given by 

vp sin <j> = — - xjj. 
And since i|/=0on the boundary, it follows that 

v n = 2irpv&m <f> 9 along the equatorial axis only, = — i// H3 



where \\t n is the value of \j* at the nth equatorial axis, or 



f^ 1 f T 2 T' ~) 



where J tt -J (Xas w ) and J' —J (X). 

VOL, OXOIT. — A. K. 
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18. The moment of angular momentum is 



hi r-nj'Z 



M = 2 2npr dr ddvp sin <]> 

JoJo 




rt f 7r/2 



a 



i//r 2 sin c/r a!0 






2tt ua H 



SiX — sin A, j 



I (J — x 2 J') x l dx 



TV .C 



St 11' 



eae 



p 



b?X — smX 



a; 



X 



sin X.t — x ,s cos Xsc — x 4 J' ) cfo 



m 



^ 



o \ 



SiX — sinX 



a 



k. 



; , W cos X 
X-y \ 



sin X s 



X 



sin X 
X 



• • • • 



(18) 



where m denotes the volume of the aggregate. 

19. The internal energy of the aggregate, supposed without translation is 






l 

9 



1 

2 



2rrp dp dz (v 2 cos 2 (j> + v 2 sin 2 <j>) 



1 \ fd^V fdf\ 2 X 2 



9 



where 
and 



7Tp 



K^J + U7 + « ; M dp<fe 



i(/ = A(J- # 2 J') sin 2 






Hence, as in the usual way. 



Si X — sin X 



E 



1 

47T 



yfr dyfr 

p dn 



I 



as 



1_ 

47T 



f * {i (7 1) 






rZ / 1 d\\r 
dz \ p dz 



vZ 1 



-y \p \ dp dz. 



Now r along the boundary \fj = 0. Also 



d ( 1 r^i/A d j 1 $i/A 



X 2 4 J sin 2 8 



cv 



P 



X 2 l|r 

d l p 



therefore 



E 



x 2 

2 7m 2 

X^A 2 

2ira 

X 2 A 2 



oirii 



X 2 A 2 

37T« 

X 2 A 2 

37Tft 

X 2 A 2 



3?m 



*ilr 2 r rfr <i# , X 2 A J' 1 , , 7 7 . 






AnrcC 



■1 frr/2 



—pj, 



2 I (J - a 2 J') 2 sin 3 dtf <&; + J' ( J - arj') x 2 sin 3 6 dQ dx 



} 



' ! ' *™ o y /«8iu^ , pnQ , \ , 2 Bin'Xa; sin 2ta 1 , 

Q3 J — o«J ( " — - £C COS AX -r* t> — £ ~ •+■ ^ COS AX HW; 

o L V A, / ' XV; 2 X?j J 



JL V 2 
5 u 



T/ /3sinX ocosX 
3 j ^__„ __ 



sin X\ , . . sin 2X 2 sin 2 X 

x + L + Tx ~ ^ 



/ 9 \ j, 2 shrX 2 sin X cos X 



iv 5 



x- 

9 
X 2 



X 



f 



J' 2 + sin 2 X I 
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or 

The energy due to translation is that due to the bodily translation of the sphere 
4- \ the same. 

The velocity of translation is 

TT a J' — \X sin X 

a X (SiX — sin X) 

Hence this part of the energy is 

3. X 4 --^3772 
2*2*3 * 

Therefore total energy is 

Q v / X \2 



7T/j/a 



(SiX - sin X) 2 
irfjira 



(Si\ - sin X) 2 



* "" x 2 j J 2 + » sin2 x + ( r J — * sin x 

2 (* " h) J/2 ~ 3X ^ Shl X + » Sin2 X ■* ' * * • (20) 



A verification is afforded by putting X = (Hill's vortex). Then 

(SiX - sin X) 2 = - s \- X G . 
Large bracket = 5TT79T9 ^°- 

which is correct. 

The preceding formulae refer to the whole aggregate. When, however, X > the 
lowest X 2 , there are more than one component, and it will be well to give the 
requisite formulae for each of these separately. Denote Xrja by y n , where r n \& the 
radius of the ^th interface from the centre. Also for shortness let S (x) denote the 
function Six — sin x. Then 

Pn S (y n ) - S (y n ^) 



fi S (X) 



(21). 



V n IT 



fin S (y n ) - S (y n ^) 



[o n -— x n J j \^/5 



J n denoting the value of J at the equatorial axis, 

M„ = ~77^ — — : — ; ( V sin y — ?/ cos y — y* ~v ) cly 
X" (SiX - sm X) J 3/n _ 1 y J J J J X 1 J J 

= x? ^)^k(?,A \ 3ynJ ( y ») ~ Syn - iJ ^-^ ~~ y2>i sin y * + 2/2 - 1 sin ^ 

5 ' X 2 

K 2 
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But 



3y - y 



9 JX 



o, 



Hence 



JLYJ-.,-^ ■■" I " '■ ■ 



1 



m/jL. 



11 



^ S. (y„) - S (y M _0 



Hyl ~ Vl-r) - ^ "/"- j J f } - yf, Pin y n + «£_, sin y^ x \ . (23). 



20. The velocity of translation is given by 



U 



fju X sin X -— 3 (sin X/X — cos X) 



3# 



X (BiX — sin X) 
ft d 



3a (SiX — sin X) d (X#) 



{J (Xx) — x 2 J (\)} x = l% 



To see how this varies with the parameter X ? refer to the graphical construction in 
fig. 2, Plate 1. The curve J and the parabola intersect in P. If A be a point on 



Fig. 3. 




the curve (fig. 3), and B on the parabola with the same abscissa near P 3 and PN be 

the perpendicular on AB, 

fju AB 



u 



oa (SiX — sin X) PN 

fji sill (a 



/3) 



da (SiX — sin X) cos a cos /3 

fi sin (a -— /3) 
3aA cos a cos /8 

Where a, /3 are the angles which the tangents to the curve and the parabola at P 
make with the axis of x, and A denotes the area of the curve OAP.MO. 

The factor --■■--.-■--- . . is always finite, except for X = 0, and positive. It is 

oa (S'bX — sin X) 

then easy to see in general how the velocity alters as the parameter X increases. 

As P (fig. 2, Plate 1), travels along the curve, U is positive. Leaving out of sight 
for the present its value for \ small, it later on diminishes to zero when P reaches a 
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certain point Q where the parabola touches the curve. It then changes sign and 
remains negative until P reaches another Q point where the parabola again touches 
the curve, and so on. 

We shall call the values of X corresponding to the Q points the X 2 values^ and 
denote them in order by X 2 1} , X 2 2) , . , . Xf/ } . Thus for values of X < X 2 J) the aggregate 
moves in the direction of the rotational flow up the axis. At X = X 2 l) the aggregate 
is at rest, the velocity of the fluid on the boundary is zero ; as X increases beyond 
this, the aggregate takes on another layer with primary rotation in the opposite 
direction, and it moves in the fluid in a direction opposed to the rotational motion of 
the innermost layer. It regredes relatively to this. The velocity at first increases 
and then diminishes until P reaches the second X 2 point, when the corresponding 
aggregate is at rest in the fluid, and so on. 

The periodic nature of the aggregates is thus evident. We get for example a 
whole periodic family of aggregates whose peculiar property is that they remain at 
rest in the fluid. The members of the family differ, amongst other things, in the 
number of independent layers each possesses. 

So we get another family formed by values of X, corresponding to points where the 
J-curve cuts the axis of x. We will call values of X, corresponding to these the X x 
parameters, and denote the orders in the same way as for the X 2 parameters. As we 
shall see shortly, the distinguishing property of this family is that in each of them 
the vortex lines and the stream lines coincide. 

For small values of X it is preferable to express the value of U in terms of the 
lowest powers of X. 

It is easy to show that 

SiX — sin X == 5-— X 3 — ... — (~~) n 



whence 



3.3! ■/ ' ' v ; 2n + 1 (2n + 1)! 



•U = 5 a (! — rfy^ 2 )- 



This gives for X = the value of U already known for Hill's vortex. 

The curve y = U/U , where U is the velocity of the non-gyrostatie aggregate of 
same cyclic constant and volume, is drawn in fig. 3, Plate 1, up to X 2 6) . The periodic 
quality is evident. 

2 1 . The directions of the lines of flow and of the vortex lines are given by 

, vp sin 6 X ijr . N 

tan © = - ~ = ~ - r . . . . . . . . . . (24 J, 

r vp cos <p a cl\jr 

dn 
X J sin 2 



cop sin y A, 

tan x = — — A = -j 

^ cop cos x a ^ 



dn 



(j - a* J') sin 2 6 
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Hence 

"X I 1 ^1 Tl ft 

tan ^ z= — tan </> — ~ '"d~T ' 1 * • (25)^ 



also 



a A J (J -^JOsin 2 ^ 

cln [ 



tan </> _ 1 2 ^ 

tan % J 



points at a less distance from, the centre than r = ~ a. At this distance x = 9 or 



Equation (25) shows that when J' = 0, i.e., for the X 2 parameters, the stream lines 
and vortex lines coincide. (It is to be remembered that we have supposed in the 
foregoing that <f> and x ^ e on opposite sides of meridian lines, and therefore 
tan <jj> = — - tan x nieans that they lie on the same side and coincide.) 

I. From X = up to X = X^, J > x*J' and J — x 2 J'< J. Hence between these 
limits, the stream lines and vortex line are on the same side of the meridians, and 
X></>> i.e., the stream lines lie between the vortex lines and meridians. At X = Xf } 
they coincide. 

II. Between X^ and X 2 1} 3>x l 3\ but J — a? 2 J'> J. For any given X, J changes 
from + to — as x passes through the value Xx = X^. For this value of x 9 or 

r = — a, x = 0« Thus, for an aggregate whose parameter X lies between the first X 2 

A. 

and X 2 roots, the vortex lines lie between the stream lines and the meridians for all 

X 
X 

the vortex lines coincide with the meridian planes, and beyond this distance up to 
the boundary the vortex lines and stream lines are on opposite sides of the meridians. 
For values of X between the first and second X 2 parameters we have to deal with 
two layers. In the outer J — • x 2 J' is negative, whilst J is negative between XP and Xf\ 
positive between Xf ] and X 2 2) . Referring to fig. 2, Plate 1, let the point jd where the 

parabola cuts the J curve be given by X', corresponding in the aggregate to a distance 

A/ 
from the centre \x = X' or r = — - a. It is clear that J (X) and J (X') are of the same 

sign. Hence, if X lies between X^ and Xf } (corresponding to P between Q 1 and R 2 ) 5 
X' lies between X^ and X 2 1} , whereas if X lies between Xf } and X ( 2 2) , X 7 lies between and 
X^L—or, taking closer limits still, between rr and Xfl We find, therefore, the follow- 
ing results. 

III. P between Q x and R 2 . In the inner spherical nucleus the vortex lines lie on 
the same side of the stream lines as the meridians— they are, in fact, exactly similar 
to the second category. At the boundary between the central nucleus and the outer 
layer <f> = 0, the stream lines coincide with the meridians, In the outer layer the 
stream lines lie on the other side of the meridian, with the vortex lines beyond. When 
P coincides with R 2 or X is the second X 1 parameter the stream lines coincide with the 
vortex lines again, but on the opposite side of the meridians. 

IV. For P between R 3 and Q 2 , we get still two layers, the boundary being given 
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by (say) X' (P atp), where JV and JX are both positive,, J — x 2 J' is positive between 
and x = X'/X and negative between x = X'/X and I. In the inner spherical nucleus 
(r = to r = X'a/X) the stream lines lie between the vortex lines and the meridians 
(similar to the first category). At the interface the stream-lines coincide with the 
meridian. In the outer layer the stream lines and vortex lines lie on opposite sides 

of the meridian for points whose distance from the centre are less than -~ a, or greater 

X 

x (2) . . \\ x (2) 

than — a. For points at a distance — a and -~ a 9 the vortex lines coincide with the 

X XX 

meridians, and between them the two lines lie on the same side of the meridian. In 
the same way the behaviour for aggregates whose parameter is greater than Xf } , may 
be determined. The periodic nature of the aggregate is again very clearly seen. 

It is perhaps easier to describe the nature of the changes above indicated by supposing 
our eyes placed in a prolongation of the polar axis. Call the vortex lines blue lines and 
the stream lines red lines, and suppose for X small that the stream or red lines lie on the 
right of the meridians. For X == 0, or Hill's vortex, the red lines lie along meridians 
and the blue lines perpendicular to these, along parallels of latitude. As X increases 
the red and blue lines swing round towards each other, the reds to the right and the 
blues to the left, and this goes on with increasing values of X up to X ( / } , when they 
coincide. Beyond X = X^ and up to X == X4 l) the red and blue lines interchange their 
relative positions. In any given aggregate the blue lines move more and more 

x (1) 

towards meridians as we pass from the centre outwards. At a distance — a from 

X 

the centre the blue lines all coincide with the meridians, both red and blue lines are 

x (1) 

swinging round to the left. Beyond the distance -—■ a the blue lines cross to the left 

of the meridians and the red lines close up towards the meridians until at the surface 
of the aggregate they coincide with them. 

Between X^ 1} and Xf we have doublets. The aggregates lying between Xi !) and X^ 
and between Xf } and Xf } are however essentially different. 

In the first set in the central nucleus the blue lines lie to the left of the red, and 
both to the right of the meridians for points near the centre. As we pass outwards 
from the centre they swing round to the left, the blue lines swing past the meridians 
whilst at the surface of the nucleus the red lines just reach it. Beyond, in the 
outer layer as we pass out, the blue and red swing further to the left, and later at 
least the red swing back again towards the meridian, coinciding with it at the 
surface. When X = \f } red and blue coincide everywhere. They lie to the right in 
the inner nucleus and to the left in the outer layer. 

Between Xf } and Xf we get aggregates in which red and blue lines again change 
sides. In the inner nucleus both lie to the right of the meridian, blue furthest out. 
They close up to the meridian as we pass out from the centre to the nucleus surface. 
In the outer layer the red lines swing further to the left and back again, the blue 
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lines follow after in the same way, crossing the meridian twice ; once in each 
direction. 

Beyond X^ 2> we get triplets* 

In general, between \f ] and \{ nJrl) the blue lines lie to the right of the red or the 
opposite according as n is even or odd. They coincide for the \ x parameter* Also, if 
n is even, both lie to the right of the meridian for the inner nucleus, the reds to the 
left for the second layer, to the right for the third, and so on* Whilst the opposite 
takes place if n is odd. 

The forms of the spirals may be obtained by finding the polar equations to their 
projections on the equatorial plane. Let (p, rj) be the polar co-ordinates of a point 
on the projection of a flow ; (p, 5) of a vortex line lying on a given sheet \jj. Then 



Pi's 



tan <j>. 



ds 



tan V| 



where ds is an element of a meridian curve* Hence 



drj 



X ^r ds 
a dyfr 

9 dn 



X -fr dr 

a dty 
9 rl0 



Provided dr is not perpendicular to ds, i.e., on the outer boundary, but then \\f = 

and 7) = 0. 

X dr X P dx 

7] = -- 



drj = 



a 2 cos 6 



vt| 



cos 6 



where x x corresponds to the inner circle of the two in which the current sheet \jj cuts 
the equatorial plane. The total angular pitch of the spiral is 



X 



'Xe> 



dx 



Xl cos 9 



B A © 9 9 



. r . (Zt>), 



where x u x 2 are the two roots of 



T , x x 9XV Xty (SiX — sin X) 7 

J (Xx) — arj X = --i^----^^-^--^ — b 9 gay. 



The above may also be written 

V = i x 



"X 



X\ 



TTfia 



J -a? 2 J' 



^ 2 J' -6 



cfe 



27). 



Equation (26) enables us easily to determine the form graphically when the 
surfaces t/* are drawn. So 



s = 77 + — J' 



■# 



{(J — X 2 J') (J — £C 2 J' — &)}~*flS 2 cfa? 



e 9 e 



(28), 



$1 



the case of a spherical boundary being excepted as before. 
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For the outside stream-lines the pitch is 

tj = 2 X -7 rf«T = X. 

For values of X, however, lying beyond X^ there are several layers in which the 
stream-lines are distinct. If x u x 2> x B . . . denote the values of x corresponding to 
the interfaces of the layers, the pitches of the stream-lines on those surfaces as we 
pass outwards are 

(x x — 0) X, (x 2 — Xi) X, (x z ■— x 2 ) X, &c. 

We have seen that on these surfaces the stream-lines coincide with the meridian. 
These parts therefore produce no part of the pitch. The twist must be supposed as 
taking place in the part of the stream-line along the polar axes. It is easy to see 
that this is so by considering current sheets near the interfaces. 

We may therefore regard the physical meaning of X to be the criterion of the total 
external pitch of the stream-lines. We will return to the consideration of the pitch, 
and the shape of these lines later. 

The total angular pitch of a stream spiral on any stream sheet \ft can easily be 
expressed in terms of the volume of the fluid inside that sheet. For 

-j -j ds XJ / p 2 m-aa 

d$ = drj — — • 77 X 



p w ] cl^r X (Six — sin X) 
dn 

j \J' 27r p ds dn . n . irua ,; 

_ ^ _ _ .__^ n ^ „ _____ y e 

Integrate round the stream surface 

XJ' cl f j 7 \J' dm. f x 

s = 7) — 7—7,-777 2rrp as dn = 77 — 7—- -777 ..... (29), 

where m denotes the volume inside xjj. 

22. The discriminating properties of the X 2 and X 2 parameters make it important to 
determine their values. The X : parameters are the roots of the equation 

J (X) ^ -7— — cos X = 0, or tan X = X 

X 

The large roots are clearly nearly (2n + 1) ~ " 

Put 

X = (2n + 1)— — = 7/ = ot — 7/ Ba y« 
Then 

COS 2/ 

— sin 2/ = 0. 
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Expanding this in powers of y, it is easily proved by successive approximation that 

1 . 2 . 13 

a oa' lb a'* 



or 



1 v ^ " ' 2 (2rc + 1) 7f 3 (2rc + 1)V 15 (2*i + 1) V 

- K . A , ft /n . 1X "63662 17201 -03558 , . 

= 1<5?079 < 2w+1 >--^ • • ■ W 

The first root is by numerical calculation 



X — 4-49341 = 257° 27' 10 



// 



The foregoing formula gives for this case (n = 1 

\= 4-49366. 

For higher values the formula is correct to five places at least 
The first three roots are 

4-49341 =s 270° — 12° 32' 50' "* 



7-72528 = 450° - 7° 22' 27" > 

- 5° 14' 23" 



a « s? a e » 



10*90408 = 630° — 

The \ 2 parameters are roots of the equation 

1 X 



cot X = 



9 > 



The large roots are clearly nearly ntr = n?r — ?/ say, where 

TlTT — ^ 1 



or 



cot y =l_ 

O 7?,vT — ■ ?/ 



A i7r — 2/ 1 

cos «/ = — r — — ) sm y« 

\ 3 nrr — *?// ^ 



(31). 



Writing nv = /8, and expanding in terms of y it is easy to prove, as in the former 
case, that 

y = ^ + * (i) 3+ * (i 

Q ' Q \3 / Q \5 

A = n?r — -§- 4 — / "■" •£ — ) 

•95493 '29026 45881 
— 3'14159n — — -■- « — —7— . , „ . . (32). 

7?y ?l 7l V ' 
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There is no root corresponding ton— 1. The first root is 

X 2 = 576346 = 360° - 29° 46' 41". 

The formula gives for this root 

\= 5-76448. 

For n > 2 it is exact to five places. 
The first three roots are 

576346 = 360° — 29° 46' 41" 

9*09506 = 540° - 18° 53' 29" 

12-32296 = 720° — 13° 56' 48" 



. (33). 



23. Equatorial Axes.— An equatorial axis is the line of particles which remains at 
rest. It is given by the equation 



or 



•by 



df 

dr 


o, 




when 


8 = 


o, 


dJ 

dx 


— 


d 
dx 


(x*J') -- 


- 0. 





The positions of the axes are, therefore, readily observed by means of the graphical 
construction in fig. 2, Plate 1 . They depend on the abscissse of points for which the 
tangents to the J curve and the parabola are parallel. For values of X > X^, the 
inclination of the parabola to the axis of x is always small. Hence the equatorial 
axes must always be near the crests (or bottoms) of the J curve, i.e., near values 
(2m + 1) J£ir. 

The equation for the axes becomes, if y be put for Xx, 

cos y+(y- jj sin y - 2f ~ ■ = ..... . (34), 

in which the roots < X are required. 

As the values of the secondary cyclic constants and other important properties 
depend on the position of the equatorial axes, it will be necessary to determine their 
values. We shall do this (1) for the case of X small, and (2) for the case of X large- 
As, however, the case of the X 2 values is special, we shall treat these separately. In 
the case of values other than X b say, e.g., X 2 parameters, all the axes of any aggregate 
depend on the particular X value. In the case of \ l9 however, they are independent 
of the particular X l9 In fact, the successive X 2 aggregates may be built up by taking 
any one and putting outside of this a suitable vortex shell. Moreover, the values of 
the axes for the X x roots are the crests, and bottoms, of the J curve, and so are 
important for their own sakes. 

L 2 
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Case of\ small— Kere y is also small If equation (34) be expanded 
y and X ? there results 



11 powers of 



X 5 



%~-2 



2\ 2 y% ( — )* 7 - - -- To - --- + J 

^ x v ; (2n + 3) (2w + 1) ! 1 

Dividing by 2y 2 /15 9 this may be written 



42,., 



u* 



?r 



r (2^"+ " 1)1 y 



YM 



0. 



X 3 



f — ~ 4 r + 302 2 (~) 



?& -f I 



7 (2k, + 3) 



i 



;(n+ i.)?/ : 



X 2 ' 11 } 



whence y can be expressed in terms of X by successive approximation, 
be found that 



y 



y 



Y 
x 



~ j 



9, 



1 



112 

X a 
224 






22 7 
5 4 



X 



9. \ 9. I 



224:, 



To A 6 it will 



(35 



This gives the equatorial axis at 



a 



^* 



v/2 L 



X^ 



X' 



v >. \ f >. 



% % 1 
224 54 1224/ J 



en X = 0, this agrees with Hill's vortex, 

'ase of \i_.~ The equation in y for this case becomes 



cos y + ( y 



1 



?/ 



sin y = 



y is always nearly n?r = t&tt — # say, where ^ is small Then 



1 



cos z — mr -— % 



sin 2 = 0, 



%tt — # 



Whence 






o 



y = ri7r 



= nrr 



nir 3 (mff 15 (?i7r) 5 
*31831 *05375 -01590 



This formula gives for the two lirst roots 

275363, 



b iiuo^j. 



The values obtained by numerical' calculation are 



274371 



6-11676, 
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The roots beyond this are therefore given by the formula correct to five places. 
The radii of the equatorial axes are r = ya/X. Hence using the values of Xj given 
in (31), the first three are. For \[ [ \ 

2-74371 



4-49341 



a = -61062 a. 



For ^ (2) 



1 5 



To 



Vx = 



6-11676 

7-72528 

2-74371 



7-72528 



a — 79179 a 



a = '35516 a 



> 



~j 



For \i»> 



9-8166:: 



o 

) 



n 



/ o 



10-90408 

611676 
To-90408 

2-74371 



r, 



10-90408 



a = '85442 a 



a = -56096 a 



a = -25162 a 



V 



Case of X large.— The number of equatorial axes depends on the order of the X 2 
parameter next greater than X. If A, lie between XSf ~ 1} and XSf , there are n such axes. 
It seems then natural to refer the magnitude of X to X£°. Suppose then 

x = xr> - x, 

where the maximum value of X is about tt,— or we may write X = X^" 1} + X, and if 
both be allowed X will have a maximum of the order ^tt. 



X^ = 



nrr 



o 

O 



^7T 



JL 

3 



3 \2 



JL 
5 



3\ 5 

7T / 



The equation in 7/ is 



cos y + ( y 



_1_ 

2/ 



sin y 



V f = o, 



in which the first n roots are to be determined. For small roots the parabola of fig. 8 
is almost coincident with the axis of x 9 and consequently the small y roots are very 
nearly equal to the corresponding values for X 1§ It will be best to obtain an 
expression for the large roots and then see how far back it holds for the smaller roots. 
Clearly y is always near mrr where m is an integer < n. 
Put 



y = rmr + z = a + z say. 



Then 



y \ yj 



v 



J (X) may be either + or — , it is of order of magnitude 1 at most. 
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Since z is not large (it is of order 1/a), we get 



1 



a 



1 



Writ 



e 



z 



a 






Zr 



cc 






i- + ,u I + ^ 



X u 



,-yo 



6 






2(-) 



VI 



a 



J'X 



X X 



a 



9,*\ 
' JL \ 

a J 



( X m O 



b 

a 



ct 



X 



z 



1 -f- ) ■":— z?z 0, 



a J X 



The greatest value of a/A is < 1, J 'A/A is of order 1/A, therefore at least of order 1/a. 
Hence in the most unfavourable cases h is < 2. The above equation can be written 



z 



z 



% 



1 

a 

b 

a 



b 
a 

a. 



bz 

a 

6- 1 



sv 

/J.-0 



£C 



3- 

9, 



a 



a 



zf 



i) 9 



33- 



6 



B***"* *" '* isuJU^n 









a;.. 



SL 

5 a 2 '* a 3 ! a.4! ' 5! 9 



(1st approx.) 



5 - 1 . (* + 2) (6 - 1) (6 - I) 2 (?) - ])* 

+ . + - s -t- • c~T~" > (^ n d approx,) 



DC 



O? 



Q 8 



«-• 



5-1 (5-1) (5 + 2) (5 + 5) 



a 



6 



oc 



6^ 



O 9 O 



e * 



. (36). 



It will be convenient to put 6—1 



v-/« JL XX Oil 



z 



c l_ c ( r3 + 3 ) (g ± 6 ) J- f c ^ + 3 K C + 6 )( c " + 4e + 6 ) 

a 6 a 8 * a' } 



o ,:i 



3e 



2 ,,3 

3 ° 



41 



c° 



5! ' 



If A is a Ai root , c 



l and 



1 



% 






7° 



15<~ 5 ? 



># 



which agrees with the result already found, 

2 4. The Spiral Forms taken by the Lines of Floiv and Vortex Filaments.— The 
equations determining these are given in § (21), Unfortunately, however, they are 
not integrable in finite forms. 

We give a graphical method for the stream-lines later. At present it is proposed 
to determine (1) the forms of the stream and vortex lines when A is small, (2) the 
pitch of the spirals near the equatorial axes, and (3) the pitch of the same on the 
outer surface. 

Let the stream surface t//, the streams and filaments on which we have to investi- 
gate, cut the equatorial plane in circles given by r/a = x Y and a? 2 . Then 



V 



X 



v; 
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where 



b = — (BiX — sin X) xh: 



In determining the vortex filaments we will take ? to be measured in the same 
direction as rj ; that is, to the right of meridians as looked at from the polar axis, 
In this case 



'X 



X 



J sm 2 6 dx 



X 



V CtOD 



Xi sin 6 (J - a? J*) 2 sin 6 cos ^ 2 } Xl y/{(J ~ of J') (J - xW - 6)} 



or 



X 



^ + 17" J / 



(1.) Case of X small. 



also b is of order X\ Put 



Hence 



X 



77 



9 



2 



ft- tvx 



</(,j - aj*j') (J _. a^J' - 5) 



J 2 T' 1 \ 4 



4,„2 



iC 



sx 2 



10 



+ 775T( aj4 "" 1 ) 



& 



30 A G. 




r 



re 2 ( 1 



X 2 



^) ( 1 - gg a? 2 + 1 



1 



V 



* s ) ™ w ( ] " a;l ) 



1- sc cfo. 



v~ 



c i 

J 



a?!, x 2 are the roots of the denominator equated to 0, viz. ? of 



X 2 



or 



a; 2 + c = -.- (x 6 - a?). 



28 



A first approximation is 



x 



1 ± y/ 1 - 4c 



rj or r 2 (say), where ^ denotes the smaller root. 



Let for a second approximation 



ar 



X — X -\- c 



X G ~ X 2 



Therefore. 



r i + £ where £ is of order X 2 e 

r\ + 2n£ — r, — £ + c = (2r, — 1) f. 



X 2 



(2r x - l)f = — (rj — n)„ 



28 



X 2 r\—o\ 
28 2^-1 



X 2 ^vv(^i + l) 
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Since 1\ -f ^2 ™ 1 » also rfc = c, 



£ = 



X 2 c r, 4- 1 



j » 



r n 



Hence the roots ar 



Vi 






Ah r\) -~ r, 



and 



<^2 



To 



X 2 c n ~j- 1 

28 " n-r, 



and the denominator becomes 

(,T 2 — xf) (Xo •— X ?J ) (1 



X 2 



28 28 



X 2 

,2\ 



"hence 



/*v>'.' 






(1 



X ) (I 



\ 



A, 






a:, 

x p 

4 J v 





(x- — xi) (xt 
1-y 



28 



t.v 



(y-ydQh-y) 



X") 



y, wnere y :-= x 



X 1 
28 



a? (6* 



cir 



In this put 



2/ 



#2+?/l ?/2 — ' 2/t 



*u 



2 



c 



308 0. 



So that 



2/ - Vi 



2 



(1 — cos 0) y 2 — J/ — 



2 



1 + cos #). 



, f 'rn 



7) 



X 

4 



X 



V 



./]. — y d0 P 



re' 

4 Jo 




1 



2/i 



(& 



1/i) sin ~ 



dd. 



X 
~2 



X 



2 



*4> 



^1 



Vi \ %/i — k 2 sin 3 <j> cUf) 





v / 1 — - y x E (£, </>) 9 where sin 2 <f> 



9 9 



and 



& 



2 _ ^Zl^l 



-2 2 



n 



r-i 



^ _l_ ri 



r. 



Arc r 2 + n + 2 
28 r 2 — r x 



28 r 



2 ' I 



8 csessasa ^[1 g* &s*s 



O 



>\*c 



9 

4aaS 



I .»— 4. C "J™ . v/ I "**"- 4 C 



X 2 c 



28 






jO 



1 -- 4 c 1 ) 



(37), 
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At the equatorial axis h = 0, on the surface k=L], Thus k increases from 
to 1 for the various current sheets in order from the axis to the/surface. The pitch 
of the helix on any sheet is 

Pitch = X \/T^y l E. 



At the surface this is X, at the axis it is 



-\ t 

X s/\ — y . — = — \f 



1 



i'l 



2? 

1T2 



° l ^ 112/ 6 



2^2 



Since 7r/(2 v /2) = I'll, the pitch at the axis is about 11 per cent, larger than on 
the surface when X is small. \ 

The corresponding quantity for the vortex filaments is given by 



* = V + 2^J' 



x 2 dx 



x x s/ (J - %-y) (J - aw - &) 9 



By what has immediately gone before 



V = 




X (Ji/X 



X 2 

(x 2 -f 1) f -\/{(l ~~ ^ 2 ) (^ 2 ~~ ^i) (^1 "" ^ 2 )} 

___/ „_______ 

X 2 X 2 \ 

1 - "oq - "28 y j V^ U 1 - y) (^ - ^i) (3/2 - 2/)} 

d0 



X 2 



X 2 



x 28 28 

15J' /•+ 



3/ /a/ 1 -y 



and 




2\Vl-yx 



J' = ^X 2 (1 — Y5-X 2 ). 



d<p 



X 2 X 2 

(1 + 3/i) - 1k> (2fe ~ 2/i) sin 2 £ f v/ (1 



28 



28 



• 2 j» \ ' 



2 sin 2 $) 



Therefore 



5 1- 



V 



10 



where 



X 2 



2 W 1 - yi i 1 - -00 (1 + 3/1) 



■* 



C?0 



(1 — n sin 2 </>) -y/ (1 — k 2 sin 2 </>) 



28 



n 



x*_ 

28 



X 



(2/2 — 2/l) = nnV(l - 4C). 



Thus 



9 = \Vl -y 1 'E(h.^) + 



28 



5-^X 2 + 



28 



2 



2X\/1 



Vi 



n (— w, &, <£). . . (38), 



At the equatorial axis n =• 0, & = ; II = tt/2 for a half turn. 
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Thus the pitch at the equatorial axis is 



5\% 



7rX 



% \2 i " /v 



and 



^ ^ 2X v 1 — 2/0 



28 



77\ 



Vo 



a I JL 



2 



X 2 
112 



Therefore 



Pitch 



nr 



X / 
\/2 



™~" ^-~~~— • I I g^Ja^a ™™~"-*~- w — ** I saaJkss* 



112/ T X 



50 A )|1 



X 2 " 

jL _A- 2™J 



7T 



5?r 



{X 2 -f 2 (5 



_3 1. \2 

JL .L <a 



— X /2 ^ 

If X = 0, the pitch is oo , as it clearly ought to be, since all the vortex filaments 
then lie along parallels. 

The Form of the Spirals near an Equatorial Axis. 

The meridian sections of a current sheet near an axis will evidently in general 
be elliptic. To find r] it is therefore necessary to determine for an ellipse the 

value of 

clr 



cos 9" 

The following general theorem enables us easily to do this, Transfer the origin 
to any point O' in the equatorial plane, at a distance c; and let the new polar 
co-ordinates of a point P be r'.ff, corresponding to r.0. Also let x.y denote the 
Cartesian co-ordinates referred to 0'* Then 



r 2 = r ■+ o u + 2(%c, 
rdr = r'cZr' + ccfe, 



dr 
cos 



?x^r 



, r cos 



r' cos & 



f rf/ 

cos 



/ + c 



CviX/ 

y' 



For the spirals near the axis the point of interest is to determine the angular pitch. 
Now clearly for a complete ellipse, whose axes are parallel and perpendicular to the 
equatorial plane, and whose centre is at 0' 



dr f 



cos 6 f 



0, 
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Further, if the axes are a, /3, respectively in the equatorial plane and perpendicular 

to it, 

x = a sin 0, g/ = y8 cos 9 

where 6 is the excentric angle of a point on it. Hence, 



ttJu 

y 



+ 2~a COS d6 
_!L /3 COS 6 



a 



IT - 



J- 



Therefore, 



/ p lit*, \ 7T0LC X 

7) (tor halt-turn) = — - . — . 



or, 



angular pitch = 



TTClC X 



(3 a 



To apply this, it is necessary to determine the form of the current sheets near the 
axis. 

Let the co-ordinates of the equatorial axis be c, o. 
The equation to a current sheet is 



J 



'Xr 



a 



9 



7 J W 



sin 2 = constant, 



or, 



9 






* 



/Xr 



a- 



= constant. 



& = J (X)/X 2 , and r, p are nearly = c. 



X?* /Xr \ 2 

Denote J ~ & [ — I by /, and suppose it expressed in terms of x 9 y co-ordinates. 



a 



a 



J 



Refer to O'. 

Then x ~ c + I? 2/ = ° + ^ where £, tj are small. Hence, if/ now denote the 
value at O', 



(o + % y 



(c + t*f + tf 



f+T^ • £+ tyV + - 1 



cox 



2 



f £ + 2&? f^- + v 2 " 



rf/ (^f a; 



rf 2 / 

**^'**' ~/i rf*™" 'VI ~ — ' — ■*■• — ■■— 
7 o 1 *-f ™3 '/ 7 -» 



9../? 1 



^V 2 , 



constant, 



dx 



dr r ? 



dy 



dr r ' 



and clfjdr = o, for c is given by this equation. 
Denote dfjdr by/ 7 , c///^ 2 by/". Then 



dy 



i 



dx 2 

dx dy 

dy 2 



r 



9 



qvi 






x- 



9 

Mil 



f" = /", since x = r to 1st order, 



xy 

7 



~ 1 JJ~t~ 1 JJ = 0> since y = to 1st order, 



1 



9 

4>« 



Z' + t/" 



.2 



0. 

M 2 
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Hence 



T 



(/+ iz€ 2 f") = constant. 



The constant Is nearly J =z f *—> a say ; then 



/ 



c 3 



-rf + yi 



or the curve Is the ellipse 



i_ jr_ 



a 



Hence the angular pitch 

Now 



ttXg 



a 




2/ 



,2/// 



/ 



// 



$ 



clr 



a 



JXr 



a 



a- 



JX 



« 2 /' 



X 2 d? 



a 




a? clif' 



3y 



f" 



„Jfx)l 



xM^ 
2 1% 2 



2 



& 



JX 

x a ~ 



2 






l 



X 3 



- > 



X/V2/ 



Jy)- 



The angular pitch of the stream-lines is therefore 



7T 




2/ 



J(y) 



2/ 8 



w 



i? 



wher 



F 



1 

2 



J(2/) 



J (y) — y' 



J(X) 



x a 



I 

9 



Now 2/ is determined by 3 {y) — y sin y + 2y 2 



2 J(X) 



0, therefore 



sm v — 2ty ~r 

2 1 "X 2 

sm y — 3y —7- 



1 

2 



y 



1 



2 « 2 



sm ?/ 

JX 



Note for X x aggregates p 



2 _ x 

■"• 2 



f 



When the value of X is fairly large we substitute for y from equation (36) 



y zzz mir + z. 
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where 



Therefore 



% 



6-1 (b - 1) (b + 2) (b + 5) 



a 



6 a 



aEEm7r , & = 2 (-yv^ 



j9< = 



1 

2 



(a + *) 2 """ 2 



( — l) m 



sm z 



(a + z) 



b 



3 



2(-) w a 2 



1 

2 



— - L 4. I 
a 2 "I" 5 



a + z b 



1 

2 



1 

a 



2 ^ 2 2*2 



z z 

— + 



a a 



2 \ / ~2 

6 



I 

2 



1 



_ 4- J 

"J™ c 



2 



w^ + r+gn)^ 



(b - If 






1 

2 



1 



3 + I 2(6- i) r "i - „. 



5-1 



l<3f 



«' 



6 - l! 2 - 6 b - 1) 



O 



> - 3 



1. 

2 



i 



2 ^ 2 6 - T 77 7 7 + 5 



7 \ ^ + 2 

6 \ a- 



1 



2 



a 



+ 2 



2 . (6 + 2) + (Ani)(* + B) 



«■ 



i 



(6 - 1) (b + 5) 



P 



2 



2^" 



1 



£2 j„ 6& _ ] 

"2a*""" 



&+2 - 



(6 - 1) (6 + 5) a 2 _ (5 - 1) (6 + 5) 



w 



or 




b+ 2 



(&- !)(& + 5)-> 



a 



i - 



W + 6& - 1 



^ 



=^v 



6 + 2 + 



(& + 2) 3 - 9& 

2a 2 



> • 



2a s 



& = 2( — ) m 



a 



JX, 



where 



and 



x = xs?*> - x 



x?> 



n7r 



7&7T 



1 

3 



3 V 



%7T 
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and 



denote nrrhyfi; then 



< 7T = g^TT say. 



a' 



! JX 



X 5 



<-Kv)'( 



1 



/3 3 



A. 



■1 



Sill { A + --g 4" -gg . 



(~r 



m 



\ 71 



\2 



1 + 



"J 



(-) 



ft. 



m 



% 



2 



therefore 



/S ( 1 



o 
O 



-X. 



3 9 

cos ( X + — + — 



s /3 



6 + 3X 2 \ f 1 / 

^ £ 2 / 1 /s V 



^v 



S 



in a, "x~ 



O 



1 7T COS X 



+ cos 



9 



sin A 



' 3 



2X 6 + 3X~ V 
1 + J + ~flT- 



cos 



2 |j + sin X 



> 



~J 



V/3 

y8 /3\/ j ' 



6 + 2 = 2 



(-)' 



m+u 



m 



% j 



2X 

1 ~T 7," ~p 

p 



9 + 6X 2 \ 

_ \ cos 



/ 9 ° "V \ 

j c*j *j u.\. \ ft "^T* 

U + ^y sin x 



$ -A- 



For very large values of X we may neglect powers of — , and then 



pitch = tt^/2 4 1 — ( — )™ +n 



m 



n 



cosX 



For the outside shell m 



n< 



pitch = 2tr sin | 



Thus in the case of the X 2 aggregates the pitch of the outer layer is very small. 
If we number the shells backward from the outside, we write n — p + 1 f° r m : 
the pitch is 

tt^/2 -{!+( — ) p [ z — £-1— ) cos X 



n 



p + 1 N2 



% 



It is seen, therefore, that there are two series of shells in aggregates of large X, one 
in which the pitches increase as we pass inwards, and an alternate series in which it 
decreases. If X lies between a \ and a X 2 parameter (X 2 > X x ), the outer series belongs 
to the first category. If X lies between a X 2 and a X, value (X 1 > X 2 ), the opposite is 
the case. In other words, if the parametral point P in fig. 2, Plate 1, lie above the 
line of abscissae, the outside layer has a very small pitch, and those of alternate shells 
increase as we go to the centre. If P lie below the opposite is the ease* 

The vortex spirals are given by 



s _ v _ _ J (X) 



lv LbJu 



x l 



{(J — % 2 J') (J — -x 2 J ; — b)Y 
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and 



Write J — x ?J J' = / 

Let x be the value of x at the axis, so that/' = when x = x . 

For points near the axis, 

x = x + f , 

/(*) =/+ af. i + \ay. ?=/+ wr- e, 



V 



XJ' f*. 



(a* 4- f) 2 ^ 



2 9 



-, <ir) W 



f-b 



\ 



2^// 



+ r 2 ^~ + f 



a 2 / 



where (\f") means the positive value of \f" . 
To the first order, 



s — ■*? 



XJ' r* (* + £f d% 



Hence for the total pitch. 






r 



^ 



2XJ'_ r*+^ 



(a$ + f) df 



(ff - r> 



Put 



2/_ 



r 



£= £i sin #„ 
_ 4\J' pi * (sjj + ff sin 2 6>) dfl 

- £fsm 2 # 



7? 



or writing 



«?(f")\ A /(_ 2/ 
V \ a*/' 



£ 5 



S— 7} 



2/ 
4XJ' 



i4» 



2/ , 2/ 



V(- 2//') Jo 
4XJ' 



*o - ^fer, + ^7 + ffi sin 2 6 
VXl^- # sin 2 6) 



ae 



xl — -^-\ F •+- -^- E 



At the axis itself F = 0. 

S — 7J : 



4X,I'^ 



7T 



V(-2#") 2 



2/ 



-p 27TA 



J/ 2 



a?f" ' -■~a^/(-2ff") 
aA /(-2//") U/;+ ° J 

av /(- 2^o w > y x 1 

2wJy 



if J{y)>y' 



J(x) 
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2it /2y 2 J\ T \ . p T 

s -'^a~a(ir" J 4 if j2/< 



v 2/ 

where y =e hx /a. 

Spirals on the bounding surface, or interface between tivo shells. This is the case 
where the transformation (Eq. 26) fails. Taking first the stream spirals 

7 X 4r ds 

dnq = - • ~^ — - 
a p ayfr 

dn 

On a spherical boundary this is zero, except for the X 2 aggregates, in which, however, 
there is no flow at all. The other part of the stream surface is the portion up the 
polar axis. Here ds = dr and dn = rdd. Therefore 

lwist on axis alone — - 



= — ar — X, 

a Jo 

There is no twist on the spherical boundary. Hence 

Angular pitch of stream spiral = X* 

Next for the vortex spirals. Here there are two portions as in the former case- 
the polar axis, and the spherical boundary. 

, X J sin 2 9 1 

P j- {(J — % 2 J') sin 2 0} 

Hence, supposing at present we are dealing with a singlet only 

2\ f a ± J 2\ r^ JX.adO 

__ ^ ^ r _j_ & _ _.. 

1 a sin 7 
o ar 



a sin 6 -J- (J — arJ ) 



r=a 



r 1 j 2\jx { irc _iL 

I T i2T' *° 7 1 * /) 

JO" ~~ 00 v Cl , y 9T/\ Jo Sill 17 

(J — • wo ) ce=:1 



Both these integrals become infinite at the poles. We must therefore treat this 
part separately. 



9 = X 



r 1 -* J 7 , 2\J(x) 



o 



j _ ^aj/ — i x gin x _ 3 j ( X ) J 



'** ^ ' X r& a 3 sin 2 6 ds 



tt sin ^9 ' a| .JA^'Y , /^V 



t . p [\ dr J + \rdOj 



> ? 
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in which £, a are small, and in the third integral, x and 9 are nearly 1, respectively. 
Let x = 1 — £ , sin = 77. Then near the pole the rectangular coordinates of a 
point referred to the pole are connected by (if/= J — x 2 J') 



/sin 2 = small constant = /3 (say), 
df 



dx 



. f sin 2 6 = /3 ? ^ 



so that when ^ =: tj each = y. 
The third integral is 



Cvflft 



J say, 



2X 



J sin 3 s/ (d? + dnf) 



r sin {{dfjdrf sin* 6 + 4 sin 2 cos 2 (//r) 2 }i 

<V (d? + d-rf) 



= 2X 



# 



(1 - f) i (#W sin 2 + 4 cos 2 -^_ . f 



JL/LhJ./ 



2XJ' 

dfl dx 



The curve is given by 



P /A + 



(,(,/rf{f 



^ 2 H- 4£ s 



df + 1 /y/ 1 



1 + (d %l drjf 



ofy 



rn 



Therefore 



& a 



r 



d?7 



»?• 



9, 



& 



2 s 



1 
2 



7 



\ 3/2 



Therefore 



e 



L 

2 



JL 



Integral 



2 XJ ' 

dfl dx 

2XJ' 

dfjdx 

2XJ' 



wmm^\w 



1 + 4rflv 6 \ 7 



1 

2 



£1 



V lip 



-f- 7 8 \ r^ 



f + 77 £ 



+ r V' 



y + 477 97 



dfjdx 

The second integral is 



7 



£log-~- + log 

£1 



7 



XJ' 



log 



7)i 



dfjdx & 7^ e 



2XJ' 

dfjdx 



pi 



2XJ' 



sin# e£/7& 



log tan 







^7T 



X<T 



dfjdx 

and is nearly = 0. Therefore 



t 1 + cos 6 

IO£f 

" 1 — cos 9 



XJ 



J' -, (1 4- cos 6) 2 

dfjdx & Q i« 2fl ' 



sin 2 



XJ' 4 

Second Integral = yjjy^ log 



and the second and third together 
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The first integral is 



Jo 



J 



J ~ X 2 J 



,2T/ 



dx = X + XJ' 



•Wi x 2 dx 



o 



,21' 



#; 2 J 



Now J — a? J' = x 2 (1 - 
does not vanish. Hence 



x 2 ) F (x), where F (as) is finite for x between and 1 and 



Int. = X + XJ' 



Wi 



dx 



X 



o (1 



-f~ XJ I 

Jo 



#) (1 + #) Fa; 
1 &; 



d[//dte 1—03 






see 



XJ' 



Therefore 



X + yrrr; log & + finite quantity, 

CvT I wX 






XJ' 4 

^ + T?rr l°g ~ + finite. 



s = )i + 



XJ' 

X sin X — 3 J' 



1 4a \/2 

log? -J- 




(1 — x) 



CLX* 



-J 



where s is the distance from the pole of the point at which the stream sheet \jj cuts 
a line joining the pole to a point on the equator. The angular pitch is therefore 
infinite at the surface owing to the filaments being parallel to the equator at points 
close to the pole. 

25. Graphical Methods.— -The graphical construction indicated in § 17 affords a 
very convenient method of obtaining a general qualitative view of the properties of 
these aggregates. It serves also for a rough quantitative one, and at least gives for 
many determinations the rough starting point which is always the most troublesome 
obstacle in numerical approximations, It may be well, therefore, here, to collect 
and enlarge on what has gone before in this respect. 

The first thing is to trace on a large scale the curve y — J (X) where X is the 
abscissa. This is very easily done, since J is expressed in simple functions which are 
tabulated. The curve is drawn for the first three undulations in fig. (2) 5 Plate (1). 
Now X determines completely the nature of the aggregate (except its volume and its 
intensity). The point P on the J curve, corresponding to X, we will call the para- 
metral point. Draw through P a parabola touching the axis at the origin. For all 
points beyond the first few undulations a circle will suffice, or the curve drawn by a 
thin lath bent to touch the axis at and to pass through P. If x denote r/a, \x 
will correspond to a point on the J curve between O and P. If P l3 P 2 denote the 
corresponding points on the J curve and the parabola, the value of i// in the aggregate 
at the point (r == xa 9 6) is given by P X P 2 sin 2 (note P 2 Pi will be negative). The 
velocity of propagation will depend on the angle at which the parabola and curve 
intersect at P (see fig. 3, § 20). If they touch, the angle is zero, and the translation 
velocity zero, In fact the parameters of the points are the X 2 values. We will call 
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them the Q points. They are easily formed by fixing a lath at O and bending it to 
touch successive loops of the J curve. It is easy to do this correct to two decimal 
places, when numerical calculation will carry it to any degree of approximation desired. 

The points where the J curve cuts the axis of x correspond to the \ x parameters. 
We will call them the R points. 

Denote the points where the parabola through P cuts the J curve again by the 
letters p. These points give the sizes of the shells into which the aggregate divides. 

If ON be the abscissa of any such point, Xr = ON, and r = -— . a gives the radius 

of the corresponding interface between two shells. It is evident at once from the 
construction that the thicknesses of the shells, as we pass in or out, are alternately 
greater and less — that there are two categories, in one of which the thickness 
increases as we pass in, and an alternate series in which it decreases. There will be, 
however, some irregularity in the two inner components. 

The position of the equatorial axes is determined by those abscissae, for which the 
tangents to the J curve and the parabola are parallel, They are easily recognized by 
the eye, and thus a starting point for calculation is readily obtained. The difference 
of ordinates of these points (PiP 2 ) is proportional to the secondary circulations of 
the corresponding shells. In fact, when multiplied by 7^/^/(S^X — sinX), the products 
give the values of those constants. It is therefore clear from the figure that these 
circulations are in opposite directions alternately, and that we get two alternate 
series of ascending and descending values. 

The function S (X) = SiX — • sin X denotes the area between the J curve and the 
axis of x up to the point X. It is clear, therefore, that it has its maximum values at 
the odd \ x points, and its minimum at the even ones. 

The tracing of the current sheets is particularly easy from the fact that they are 

given by functions of the form 

xfj = F (r) . sin 2 0. 
Let 

a i at 



/( r) =j(^\_jij ( x), 



and let \jj and r denote values at the equatorial axis (i.e., i/r a numerical maximum). 
Then 

?Y~ := "77 r Sill u, sin U °zz A / — — , — ---, 
fo f(r ) V f f(r) 

On squared paper, draw a series of circles, radii sub-multiples of a, say at intervals 
of *05a or 'la, also the circle r = r . This last circle has the property that all the 

X X «/ 

current sheets cut it at right angles. 

Let us trace first one sheet (say xjj = 'l\jt ). We do this by tabulating the values 
of sin 6 for values of r, corresponding to the series of circles drawn. Now mark on 
the bounding circle (r = a) points whose abscissae are those tabulated values (which 

N 2 
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is done at once on the squared paper). Mark the points where the radii vectores to 
these points cut the corresponding circles. Join these points by a continuous curve, 
and the shape of the particular xjj curve is obtained ; call it the \jj x sheet. This first 
curve should be obtained with care and as much accuracy as possible. We may now 
proceed to draw from this as many of the other sheets as we please. Suppose we 
want to draw the curve \jj = k. \jj . We set a pair of proportional compasses (or any 

similar method) to the ratio \/lOk. Suppose the \p t cuts any particular circle at P, 
set the short legs of the compasses to its abscissa-. Turn it round and find the point 
on the same circle whose abscissa is the new value. Proceed thus with the other 
circles and the sheet is rapidly traced. Although this may appear cumbrous in 
stating, it is very expeditious in practice, and with a moderate amount of care very 
accurate. 

Having traced the \js curves, we may now easily trace the proj ections of the stream 

lines, for these are given by 

fcLx 
^~ Q = 2PN (see fig. 4). 

Fig. 4. 




\ 



\ 



26. It will be interesting to go into further details for a few cases, and for this 
purpose we take the first two aggregates of the X 2 and X x families, 

The distinguishing feature of the X 2 types, is that the aggregates are at rest in the 
surrounding fluid. The distinguishing feature of the X x types is that the vortex and 
stream lines are coincident. 

X 2 aggregates. Here 

mfjb ~~ X sin X 



U= 0. M 



E 



« ^ A 



15 Six 
X 2 sin 2 X 



sinX 



(SiX — sin X) 2 



7. 
9 



X 2 sin 2 X -pp 



where E is the energy of a Hill's ,i o y/egate of equal volume and circulation. 



PROFESSOR W. M. HIOKS ON VORTEX MOTION. 



93 



The first parameter is X^ = 5 "7637 = 330° 14'. 
The equatorial axis has a radius = *5130a, 

M = '09$5mp 

E = 1-6979E 

v— 2-1020/*. 

Angular pitch of stream lines at surface 

j j j 3 at) clX IS 

„ vortex lines at axis 



5? 

5> 



330° 14'. 
334° 58'. 

267°. 



The forms of the current sheets (i//) are drawn in fig. 2, Plate 2. The projections 
of the stream lines in fig. 3, Plate 2, These latter were determined by the graphical 
method described above. 

The second X 2 parameter is 

\f = 9*0950 = 3?r — 18° 53' 40". 

The equatorial axes are given by 

\ 2 x = 2*6616 and 6'2718, 
or r — *2926a and # 6896a. 



Radius of internal nucleus = *4694a, 



M = 

Tf 1 — 



*1459m/x 

3-727E 

1-9403/a 

•9403/x 

2'063 



M 2 



Vo 



•08904m,/x 1 

*1890m 2 /x 2 

1'1747/xi 
3*641 5jLt 2 . 



5? 



?? 



55 



5? 



35 



Total angular pitch of stream lines outside 

inside nucleus 
outer shell 
Angular pitch of stream lines at 1st axis 

„ 2nd axis 

vortex lines at 1st axis 

2nd axis 



5> 



?? 



?? 



? J 



521° C'. 

244° 37'. 
276° 29'. 
284° 21'. 
320° 9'. 
308° 48'. 

A: a A I L . 



The Xi aggregates. 

U 

M 



JL.I.fcJI. t5 



V 



sin X 



3a BiX — sin X 
— sin A, 






W/X 



X (Six — sin X) ' 
sin 2 X 



sin 2 X 



— iTu a — ^ — *— - —™ ■ — ; — 

9 ^ (S*X-sin\)*~ 9 (S^X-sinX) 



E 0! 
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where U is the velocity of translation of a Hill's aggregate of the same volume and 
circulation. 

The stream pitch of these aggregates at the axis takes a very simple form, viz., 

V = 



aJ[ 






The first \, root is \f = 4-4935 = 257° 27' 30." 
The equatorial axis has a radius — "6106a. 

U = -6189U0, M = -0826m/*, 
E = 10'724E , v— '8016/x,. 

Angular pitch of stream lines at surface = 257° 27' 30". 

axis =297° 4'. 



rO At 
IA Y I m — — v. *-* / 

5? ?? 35 



The forms for the stream sheets i// are shown in fig. 4, Plate 2. It is to be noticed 
that there is a considerable difference between the angular pitches outside and on the 
axis, whereas in the X^ aggregate they were very nearly the same. 

The second \ x parameter is Xf> = 77253 = 450° - 7° 22' 27" 

The equatorial axes are given by 

y = X lX — 2-7437 and 61168, 
or, 

r^r.'S552a and '7918a, 

Badius of internal nucleus = '5816a. 

U = — 3'094U 
M = — -2403w/x, 
E = 26'803E 
fjLj. = 4'9740/x, 
jx 2 = — 3'9740/x, 

Total angular pitch of outside = 442° 37' 33". 

„ „ on inside nucleus = 257° 27' 30". 

Angular pitch at inner axis =297° 4\ 

I, ,, oncer 3 , ^^ Zui £>y » 

In all the \ x aggregates the expression for the angular pitch at an axis is 



8 S it9 tHMHMMCW 


— 1-2500 


Mi = 


•0826m,/t 1 


M 2 = 


•1005m 2 /i 2 


<■» a eBSKsBssoa 


l-2707/i! 


I> 2 r= 


JL O .1 oO LJLq 



V(* 



7T 



y 
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Hence, when X x is large, the outer layers have their pitches at the ax:es about 

7jV2 = 254° 31'. 

Fig. 5, Plate 2, shows the relative positions of the shells and axes for the XI and X? 
aggregates. The thin lines belong to the X 2 , the dotted to the X x . A, A are the 
position of the X 2 equatorial axes. B, B those of the X^ 

27. In the preceding investigation we find doublets, triplets, &c, naturally arising. 
We may have also built-up systems consisting of monads, dyads, &c, as in the cases 
developed in the previous section. Each element of a poly-ad may consist again of 
singlets, doublets, &c. I do not propose now to develop this theory of multiple 
combination to any length, but merely to draw attention to it, and to determine the 
necessary conditions for the case of a dyad only. 

Referring to § 15, the general solution of the differential equation contains not only 

COS V 

J functions, but also the functions Y 2 = + sin y 9 which are suitable only for space 

y 

not containing the origin. They are therefore suitable for any shell embracing an 
interior aggregate. In the shell the functions will be of the form AJ + BY, or as it 

may be written 

sin (* + y) / • \ 

— — ~ - COS (a + y). 

It will be convenient to denote this by f(a, y). 

Let now the radius of the interior aggregate be a, that of the exterior b. Let 
also X, X' denote the corresponding parameters. 

Then we may write 

Inside fa = L 



- ) - ~ JX, 

a / a? 



sin 2 ...... (39), 



\V\ r 



Shell i> 2 = qL\f(u,'f)-^f(a,\')\sixrd . . . (40), 



Outside xfj — — - 7r V ( r 2 — — j sin 2 6. 



At the interface fa == fa and fa = 0, therefore 

/ \'a\ a 2 



/{"'"Tj-Vf^'^ - 



Write ajb =p. This equation, when developed, gives 

tana Y(X»-^Y(V) < 41 )' 

Moreover, the tangential velocities must be the same. Hence, when r — a, 

dxpt/dr = d\\i 2 /dr. 
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Therefore 

~ {X sin X - 3J (X)} = ^ {-/(a, X» + X'p sin (a + X'p) - Iff {a, X')}. 
But 

f(ct,\'p)=p 2 f(a 9 Xy 

Hence 

X sin X — SJ (X) = g {\'<p sin (a + Xp) — 3/(a, X^>)} . . . (42). 

So, also, making d\jj 2 /dr = d^/dr when r = 6 we get 

V = j|{/K)-*X'BinX'} (43). 

Equation (41) determines a ; Equation (42) gives a relation between X, X', jp, and g , s 
We can therefore Impress in general three further conditions. For instance, ratio of 
volumes, ratio of primary circulations, and ratio of secondary circulations. 

There Is a natural connection of the various singlets which go to make up an 
aggregate of the kind first discussed. At any interface all the differential co-efficients 
are continuous. In the polyad aggregates this Is not so. Differential co-efficients 
beyond the first are not continuous. Monads, &c e| which go to form them, are arti- 
ficially combined. It is possible we may, on this basis, develop a theory of special 
aggregates which will unite with one another, or split up and be capable of uniting 
again In another manner. Some progress has been made with such a theory, but 
before an attempt Is made to carry such a theory out It will be necessary to investi- 
gate the stability of the various systems. I hope soon to be able to take up this 
question* 

[May 6, 1898,— By the permission of Professor Hill, to whose careful reading of 
the MS. I owe a great debt, I append an independent and very suggestive proof by 
him of the general theorem of gyrostatic vortices, based on the equations of motion. 

Take as co-ordinates r 9 0, z. 

x = r cos 

y = r sin 

Let p be the pressure, 
p the density, 
V the potential of the impressed forces. 

Let r be the velocity increasing r, 
cr be the velocity increasing 0, 
iv be the velocity increasing z. 

Then the equations of motions are 
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dt dr r dO dz J r dr \ p 

d d d d\ ar d I p , 

tf£ a?* r #0 dfe / r r dO \ p 

/ d d d d \ dip 

{7u + r ^ + " Vie + w * j w = ~ & 1.7 + 

d . x da d , x 

7r ^ + 5* + * ^ = °" 

It is desired to find a solution in which all the quantities are independent of 6. 
Therefore 

dr da div ^ d f p _„.\ 

do ? d,e v > de ~ u? de Vp ^ j ~ 



The last gives 



d d d \ . ■. 

— + 7 — -j- w — (rcr) = 0. 

dt dr dz / x 7 



If therefore i// be the equation of a surface always containing the same particles of 

fluid, it is possible to take 

1 
T(r _ j ^y 

-txiSO 

i (rr) + ~ (rw) = 0. 

Let /c be the current function (which I distinguish throughout from \jj). 
Therefore 



1 d/c 1 <i/<; 



Substituting in 



it follows that 



2irr dz 5 " 2irr d? 



I d , cZ , rf \ f 



d-\\r 1 rZ/c $>|r 1 dfc d\fr __ 

dt 2irr dz dr 2irr dr dz 



Now make the further supposition that the surfaces xjj = const, move without 
alteration parallel to the axis of z with velocity Z. 



Therefore 



d^r • r/^ 

— - /a 



dt ~ cfo 

Therefore 

d4r d , • x <^ilr <i , • 

_l /^ _ nr-z) —. -J- — (k — ■ 7rr Z) 

VOL. CXCIL — A. O 



0, 



98 



PROFESSOR W. M. HICKS OIS T VORTEX MOTION 



Hence we can take 



Therefore 



We have now 



Now 



Therefore 



K 



7TV 2 L 



\fJ. 



T 



1 df 

2rrr dz 



W ~ 



1 dyfr 



9 



rrrO ' 



ITi 



dr 






d 



V 



dr \ p 



d 

dz 



V 



+ y + 



-|- v + 



9 . 9"» 



*T + ttf" 



T \ dz 



fdr 

' dz 



d: 



dt 
dw 



z dr r 

dz 



z 



dw 
: "dt 

dw 
dr 



f dr 



dz 



div 
dr 

dw 
dr j 

diu\ 

dr 



dt 



= -z 



C?7t' 



rf; 



O" 



— ^ + v + — 



\ 



Z^) = (iy- Z) 



'd 



T 



1 fltyr /^T 

2irr dr \ dz 



dur 
d? 



/ 



dw\ 
dr J 



4-Trh 



cs LfWJ 



i 



4tT 



s [/wo: 



d /> + y + T ' 2 + v} 



dz 



2 



Hence 



Therefore 



d 



fdr 



d^r I d: 
-.dr \ 



dio K 
dr 



'dr 
dyfr d [ dz 

dr dz 




IV 



fdr 

1 dyjr fdr 

2irr dz \dz 



dw 
dr 



dw N 



EM 

2rrr A . 









d dz 



dr 



dr 



dr 



I dr 
dyfr J dz 



dz \ 



dw\ ~" 



dr 



r 




fW/'W 



irr 



0, 



Therefore 




d^r d 


i 


fdr 


dr dz 


r 


\dz 


Therefo 


re 





dw\ /(*)/(*)' 



dr 



L( dT 

r \dz 



2irr 2 



dyfr d 

dz dr 



1 



dw 

dr 



'dr 

dz 



/(+)£» > + F Wi 



;7ry" 



rfw\ f {yfr) f (yfr)' 



dr I 



lift* 



= n 



dr 2 



1 cZ^ cP-v/r __ 



r 



2wr»rti,)-fWfW. 



PROFESSOR W, M. HICKS OK VORTEX MOTION. 



99 



Therefore 



d I P _i_ v j_ T " + w% 



and 



Kf+^ + 'T'-' 



Therefore 



^ 



+ V + 



9 i '> 



Ziv 



7a%v 



1 <ty 



27rr $r 



VM + ,. FW) 



-J 



clr 



8ttV* "*" 2tt 



7( ^W'« 



1 dyjr 

2irr dz 

r.fe 8ttV 8 ^ 2tt 






9 



Zw + n vr H o — = arbitrary function of i. 

87rV 27T ^ 



Therefore 



P 



1 



- -[. V + -| [t 2 + or 2 + (to — Z) 2 ] + ^~~ F- (i//) = arbitrary function of £. 



This arbitrary function of t is in this paper always a constant. 

The last equation, together with the following, are the important equations : 



* 

k = xjj + 7rZr 2 ? 

1 rZ/c 
27rr ^ 

1 die 



- = i; /WO 



W 



5 



1 <&Jr 

27rr ffe 



i^ 



1 rty 



2tt?' clr 



2irr ch 



~ ~p /j 



clr dw 
dz dr 



9 



1T7 



fty)f'W+rYW), 



d; l f 

dr 2 



f ? + SJ 1 = -/<*)/» - 2^F(# 



Whenever the conditions for the continuity of the r and w components of the 
velocity have been satisfied at a separating surface whose equation is \]j = const., 
then if the irrotational motion outside the surface have cr = 0, we must have <r = 
when i// is equal to the parameter of separating surface, if there is to be no slip there. 

Therefore f(\ji) = 0, when i// is equal to the parameter of separating surface. 

This is the case in the Third Section of the Paper." 

o 2 
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